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Abstract The first part of the paper reviews applications of 2-spinor methods to rela- 
tivistic qubits (analogies between tetrads in Minkowski space and 2-qubit states, qubits 
defined by means of null directions and their role for elimination of the Peres-Scudo- 
Terno phenomenon, advantages and disadvantages of relativistic polarization operators 
defined by the Pauli-Lubanski vector, manifestly covariant approach to unitary rep- 
resentations of inhomogeneous SL(2,C)). The second part deals with electromagnetic 
fields quantized by means of harmonic oscillator Lie algebras (not necessarily taken 
in irreducible representations). As opposed to non-relativistic singlets one has to dis- 
tinguish between maximally symmetric and EPR states. The distinction is one of the 
sources of 'strange' relativistic properties of EPR correlations. As an example, EPR 
averages are explicitly computed for linear polarizations in states that are antisymmet- 
ric in both helicities and momenta. The result takes the familiar form ±pcos2(a — /3) 
independently of the choice of representation of harmonic oscillator algebra. Parameter 
p is determined by spectral properties of detectors and the choice of EPR state, but 
is unrelated to detector efficiencies. Brief analysis of entanglement with vacuum and 
vacuum violation of Bell's inequality is given. The effects are related to inequivalent 
notions of vacuum states. Technical appendices discuss details of the representation I 
employ in field quantization. In particular, M-shaped delta-sequences are used to define 
Dirac deltas regular at zero. 

Keywords Two-spinors • SL(2,C) ■ qubits • quantization 
PACS 03.65.Ud • 42.50.Xa • 03.70.+k 



M. Czachor 

Katedra Fizyki Teoretycznej i Informatyki Kwantowej 

Politechnika Gdafiska, 80-952 Gdansk, Poland, 

Krajowe Centrum Informatyki Kwantowej, 81-824 Sopot, Poland 

and 

Centrum Leo Apostel (CLEA) 

Vrije Univcrsitcit Brussel, 1050 Brussels, Belgium 

E-mail: mczachor@pg.gda.pl 



2 



1 Introduction 

Eiristein-Podolsky-Rosen (EPR) paradox [T] and Bell's theorem [2] — cornerstones of 
quantum information and cryptography — deal with relativistic issues of locality, but 
methodology of proofs is non-relativistic. In the mid-1980s some authors, including 
myself, tried to amend the inconsistency and reformulated Bell's theorem in relativis- 
tic formalisms (relativistic quantum mechanics 3 , algebraic field theory [4]). First 
relativistic results were not entirely trivial (violation of Bell's inequality by vacuum 
fluctuations, non-invariance of EPR correlations for spins) but did not attract much 
attention. The turning point for the emerging field of relativistic quantum information 
was the paper by Peres, Scudo and Terno [5] on non-invariance of entropy defined 
in terms of 'reduced spin density matrices'. The line of research following from [5] 
culminated in the review paper [6j and is still continued by various authors. 

The goals of the present paper are, in a sense, complementary to those of [5J. I want 
to concentrate on manifestly covariant 2-spinor approach to storage and communication 
qubits, technically and philosophically inspired by the works of Penrose [3|8], but not 
widely known and not included in [BJ . One of the unusual features of qubits formulated 
in a 2-spinor way, and going beyond the standard helicity formalism, is the fundamental 
role of projections of spin on null directions in space-time, even if massive particles are 
concerned. In particular, it is always possible to project the Pauli-Lubanski vector on 
principal null directions of SL(2,C) transformations. Decoherence of the Peres-Scudo- 
Terno type is then eliminated, at least for massive particles [9]. 

The second class of results comes from application of null directions to electromag- 
netic fields. It can be shown that a class of twistor-type spin-frames has covariance prop- 
erties leading to 4-potentials that do not change gauge after Lorentz transformations 
(i.e. A a (x) is a world- vector field). Simultaneously, the same twistor-like transforma- 
tion implies that momentum-space amplitudes of A a (x) split into classes belonging to 
different representations of the Poincare group: Spin-1 zero-mass representation which 
we regard as photons, and two additional scalar massless fields. Electromagnetic qubits 
are associated only with the spin-1 part. 

The construction can be performed for both first- and second-quantized fields. 
We concentrate on quantized Maxwell fields, but quantization procedure is not the 
usual one. I believe the procedure I employ, based on reducible representations of 
harmonic-oscillator Lie algebras, is in many respects superior to the standard one, 
where irreducibility is implicitly assumed. Readers interested in more details of the 
relativistic formalism based on reducible representations are referred to [10] , 

In order to analyze relativistic EPR correlations of photons we have to make the 
notion of an EPR state more precise. The EPR paradox is based, technically speaking, 
on maximal entanglement in at least two complementary bases — I term this property 
the EPR condition. In non-relativistic quantum mechanics EPR states are simultane- 
ously maximally entangled in all bases, and maximally symmetric. In the relativistic 
case we have a similar object, the two- index antisymmetric spinor eab- The problem 
is that eab is invariant with respect to spinor representations, while these are the uni- 
tary representations we have to work with. Scalar states were used in the EPR context 
by Caban QT|, but it turns out that typical scalar states cannot be identified with 
EPR states, unless the two photons have the same momenta. The EPR condition is 
not SL(2,C) invariant if the two photons have different momenta, a fact responsible 
for relativistic non-invariance of EPR correlations for linear polarizations (and certain 
interferometric qubits) . 
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Another conceptual difficulty with relativistic qubits is the choice of yes-no observ- 
ables. For massive particles a natural definition is given by the Pauli-Lubanski vector 
(cf. [T2]) but the resulting observable is not linear in number of particles. I show how 
to define an analogue of the Pauli-Lubanski vector which is linear in number of par- 
ticles, but this is done in representations of harmonic oscillator Lie algebras that lead 
to a well behaved vacuum part of 4-momentum (this is one of the reasons why the re- 
ducible representations are here useful). Still, components of the Pauli-Lubanski vector 
commute in the massless case, so are useless for EPR correlations. 

Since the field is quantized in a nonstandard way, one might ask to what extent the 
EPR averages for linear polarizations and interferometric observables are sensitive to 
modifications of field quantization paradigm. The answer is negative: In both reducible 
and irreducible representations of oscillator algebras the EPR averages possess the 
well known cos2(q — /3) term, multiplied by p parameterized by sets of wave vectors 
analyzed by photon detectors. 

The paper is organized as follows. In Sec. 2 I recall basic links between 2-spinors 
and null vectors, and stress certain similarity between tetrads in Minkowski space and 
2-qubit states. This section simultaneously introduces abstract-index notation needed 
later in the paper. Sec. 3 deals with first-quantized Dirac fields. I analyze properties 
of the Pauli-Lubanski vector and show why it pays to work at a level more general 
than spin defined via helicity. Sec. 3.5 discusses an unorthodox approach to unitary 
representations of the Poincare group, regarded as passive versions of active spinor 
representations. In Sec. 3.7 I concentrate on the Peres- Scudo-Terno problem and show 
that a formalism based on principal null directions plays in this context a role of 
error correction. Sec. 4 introduces basic 2-spinor concepts needed for massless fields. 

1 emphasize the special role played by twistor-like spin-frames and, in Sec. 5, apply 
them to first-quantized electromagnetic qubits. In Sec. 6 I perform field quantization 
in terms of reducible representations of harmonic oscillator Lie algebras, and define 
polarization operators that are linear in numbers of particles. In Sec. 7 I concentrate 
on relativistic aspects of EPR states. EPR correlations are computed in Sec. 8. The 
main result shows that the reducible representations I advertise predict essentially 
the same formula for EPR averages as the usual formalism. This is perhaps the most 
important new result of the paper. Sec. 9 discusses the issues of 'entanglement with 
vacuum' and 'vacuum violations' of Bell's inequality. I illustrate on simple examples the 
ideas behind these approaches, and show that they are related to different meanings 
of the notion of 'vacuum'. I end the paper with technical Appendices related to the 
choice of representations employed in my approach to field quantization. 

2 Two-spinors: Geometric and algebraic preliminaries 

Let us consider three components (x,y,z) of a unit vector, 

x 2 + y 2 + z 2 = l, (1) 

taken in some basis in R 3 . ([T} describes a unit sphere whose points can be equivalently 
parametrized by means of stereographic projection on the plane z = 0. A particularly 
useful form of this map is obtained if the plane is parametrized by complex coordinates 
linked to the point on the sphere by £ = (x + iy)/(l — z). Let us note that the 'north 
pole' (0, 0, 1) corresponds to £ = oo, whereas the 'south pole' (0, 0, —1) implies £ = 0. 
Yet another possible parametrization of the sphere is in terms of projective coordinates 
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(£, rf) satisfying ( = (,/r). Now the poles correspond to (£, 0) and (0, rf). The coordinates 
are projective since (£,77) and (X^,Xr/), for any A 7^ 0, define the same point on the 
sphere. 

Let us now treat our sphere as a light cone in Minkowski space, 

1 - x 2 - y 2 - z 2 = 0. (2) 

The four components (l,x,y,z) of a null future-pointing world-vector are linked to 
projective coordinates by 

(i, x , y , z) = ( i, ii+M, -iii^zi, . ( 3 ) 

Multiplying both sides by (£f + rjfj) / y/2 we obtain coordinates 

(t, x, y,z) = ±= Ui + rfij, & + & -m - v o, ei - W) w 



V2 

of a point belonging to the same null line in Minkowski space as the point corresponding 
to (l,x,y,z). © is future-pointing since T > for a non-vanishing There is, of 

course, nothing fundamentally special about l/y/2, but this choice is convenient, as we 
shall see later. Given a null future-pointing world-vector with components (T, X, Y, Z), 
and solving the four equations @, one can determine £ and rj up to a common phase 
factor. (T, X,Y, Z) is termed the flagpole of (£,»?). A flag plane can be defined by a 
vector tangent at (x,y,x) to the sphere |T]). The flagpole and the flag plane determine 
(£,77) up to a sign. This sign ambiguity is related to topological properties of the 
rotation roup. 

I have purposefully stressed that (T, X, Y, Z) are components of a Minkowski-space 
world- vector evaluated with respect to some basis. The pair (£, rj) can be associated 
with any basis in a 2-dimensional complex space, but there is no obvious link between 
the Minkowski space and this 2-dimensional linear space. Putting it differently, if t, x, 
y, z, define a Minkowski tetrad the vector Xx + Yy + Zz + Tt can be linked with a 
vector (a 2-spinor) £o + rjL, where o and l are basis vectors in C 2 , but there is no a 
priori link between t, x, y, z and o, i. However, let us consider the tensor product of 
this spinor with its complex conjugate, 

(£o + r]i) <g> (£o + fjt) = ££o ® o + £770 <g> I + -q^L <g> o + rjfjb ® t (5) 

= (V 7 + rf) * -<4= (ft - tf) y + i (tf - z 



X 



where 



+ — ^ + v fj)t, (6) 



-i= (o ® o + t ® , (7) 
-i=(o® i + t.® o), (8) 



(9) 



4=fo®o-t®fc), (10) 
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and 



o 




(13) 



(12) 



(11) 



(14) 



A change of basis {o, l} — ¥ {o' , l'} implies an associated change of components (£, r/) — > 
(£ , r/). If (£, r;) — >■ (£ , ?/) is a linear transformation with unit determinant (an element 
of the group SL(2,C)) then the associated transformation (X, Y, Z, T) -> (X\ Y' , Z' ,T') 
turns out to be a proper isochronous Lorentz transformation. Since, Xx + Yy + Zz + 
Tt = X'x' + Y'y' + Z'z' + T'i' , the map {x,y,z,t} -> {x 1 , y' , z' , t'} is a proper 
isochronous Lorentz transformation as well. Tensor products of 2-spinors with their 
complex conjugates transform under SL(2,C) transformations in the same way as 
world-vectors in Minkowski space under boosts and rotations. This leads us to the 
abstract-index formalism. 

2.1 Abstract-index formalism of Penrose 

Second-rank spinors {x , y, z, t} have transformation properties of a Minkowski tetrad. 
Spinor approach to space-time geometry is based on the assumption that one can 
identify the two bases, {x, y, z, t} in C 2 C8> C 2 and {x, y, z, t} in Minkowski space, and 
simply skip the tildas in (I7)l- pup . From this perspective 2-spinors are more fundamental 
than world-vectors in Minkowski space. The Minkowski space can be regarded as a 
structure derived from a more fundamental spinorial level. 

This is the starting point of the abstract-index formalism of Penrose [7], where 
f7))- (fTU|) would be written as 



The indices are just labels (analogous to 'Alice' and 'Bob'; think of t alicc = t Mice < mce ') 
and do not take numerical values. The null tetrad corresponding to (|ll[l - (|14p is ex- 
pressed in the abstract-index form as 




(15) 



(16) 



(17) 



(18) 



■a 




(19) 
(20) 



m 



a 



A -A 
i L 



A -A' 
i o 



(21) 
(22) 



n 



a 



A -A' 

I L 
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In this approach there are no operations that would involve (anti-)symmetrizations of 
primed and unprimed indices. In consequence, we do not loose generality by assuming 
that the primed labels occur to the right of the unprimed ones. In particular, we identify 
a APA' with fij^/otA for all 2-spinors and Pa 1 , and one does not need to distinguish 
between hermiticity and reality. As an illustration, let us consider reality properties of 
the null tetrad: 

= = o A ' o A = o A d A ' =l a , (23) 

= &tW = o A ' i A = l A o A ' =m a , (24) 

and so on. This formalism is very convenient and does not lead to ambiguities in 
practical computations but may seem somewhat counterintuitive, especially to those 
who are accustomed to the more traditional spinor notation employed in relativistic 
quantum mechanics. This is why in the next subsection we shall describe a variant 
of the Penrose formalism that can be directly translated into formulas we know from 
quantum mechanics textbooks. 

But before we do that let us discuss relations between abstract and numerical 
indices. Numerical indices, and 1, are sometimes needed but then we denote them 
by upright boldface fonts. Accordingly, the symbol 4> A denotes a 2-spinor ('a spinor 
of Alice'), but 4> A may equal <f>° or cf) 1 . cf> A is basis independent, but (f> A implicitly 
depends on a basis. Now consider two 2-spinors, a A — cPo A + a 1 l A — a°o A + ot 1 l A , 
f3 A = f3°o A + p x i A = /3°5 A + P X l A , where the two bases are related by o A = (So) A , 
l A = (Si) A , detS = 1. The determinant 

= a p - a P = c AB a P = a B /3 = -a P A (25) 

is independent of S and defines 

,ab=(_° 1 ;). (26) 

The last two equalities in (|25p show how to lower spinor indices in SL(2,C)-invariant 
manner: = « A £ab- The inverse rule is £ AB q:b = a A , with e AB = £ab- 

The lower-index ag may be regarded as a component of a spinor a a dual to a A , the 
duality being given by oaP A = cha/3 A - The formulas ag = a A SAB and e AB ag — a A 
define at the abstract-index level the isomorphisms between modules of upper-index 
2-spinor fields and their duals. Note that e CB eca — £ A B ac ts by £a B 4>B — 4>A> and 

«■-(;;)• ^ 

ea B thus plays a role of spinorial Kronecker delta; at the abstract-index level ea B is 
the isomorphism between 'the spinor of Bob' and 'the spinor of Alice', a map quite 
similar to the teleportation protocol [13] . 

The basis o A and l A , satisfying oa^ = 1, is termed the spin-frame. Since (Pa4' A = 
for any cf>A, spin-frames possess a kind of gauge freedom: (oa + Xla)^ = °A{ iA + 
HO A ) = oa^ A = 1. It is interesting that in electrodynamics this ambiguity indeed man- 
ifests itself in a form of gauge transformation associated with Lorentz transformations 
of four-potentials. Any spin-frame satisfies 

SAB = °A^B - L A°B) (28) 
£A B = °A B - L A° B • (29) 



a 


a 1 




5° a 1 


P° 


P 1 




P° P 1 
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If g ab is the Minkowski-space metric tensor then g a b must be the abstract-index form of 
Kronecker's delta in Minkowski space. Accordingly, g a b = sa B ea /B ( £ A B — e A' B )• 
Lowering the indices we obtain the following three useful abstract-index forms of the 
metric: 

g a b = £AB£A'B> (30) 
= t a t b - x a x b - y a y b - z a z b (31) 
= n a l b + l a n b — fh a m b — m a m b . (32) 

The reader may have noticed that in the Penrose formalism the Minkowski tetrad 
has all the properties of the Bell basis of two-qubit entangled states. The null basis 
I — o®6 = l a ,m — o®l = m a , m = l®6 — fh a , n — l®Z = n a is analogous to the basis 
of product states. There exists an unexplored 'duality' between space-time geometry 
and quantum information theory. Some preliminary considerations in this spirit can 
be found in [13], where analogies between teleportation protocols and metric tensors 
of Lorentzian manifolds were discussed. Lorentzian techniques, with applications to 
multi-qubit entanglement, can be found also in |14H15|[T51ll7| . 



2.2 Abstract-index analogues of Pauli matrices: Infeld-van der Waerden tensors 

Let us return to the formula Q, but written as 
X = (A , X 1 , X 2 , X 3 ) 

= -t=^°0° + cp 1 ^ 1 ^"(p 1 + ^(p , -i^ ^ 1 - (^(jP ),<p°<j>° - (p 1 ^ 1 ) (33) 
I 1 2 3 u AtA' , Q/I n 

= (g AA',g AA',g aa',9 aa>)<p 4> , (34) 

A a = } a AA^ A ^ A • (35) 

The coefficients <? a AA' are known as Infeld-van der Waerden symbols, and can be 
grouped into four matrices 

o _ (g°o > 9°0l>\ _ 1 f l o\ A A' r „ fi s 

,1 / g 1 00' 9 01' \ 1 (0 1 \ „ AA' 



" AA ' = 1 -\o>g\v) = T2V*) =gi ' (37) 
I g 2 oo> g 2 oi> \ _ i / o -A _ a a' 



z _ i a uu' a ur i _ _j- i " • _ „ /qo\ 



3 3 00' ffV \ = _L f 1 A = AA' 



„ 3 I y uu' a ui' \ * i * " \ a a / Qn \ 

9 AA ' = {g\o>g\r)=7t{°-i) =93 ■ m 

The numerical indices are lowered or raised by means of the epsilons and t/ ab = g a \j — 
diag(l, — 1, — 1, —1). In quantum mechanics textbooks one works with matrix 4-vectors 
cr = (l,a x ,uy,a z ) and a = (1, — a x , — — Oz) = eae T . Infeld-van der Waerden sym- 
bols, with appropriately raised or lowered indices, play precisely the same role, but 
make the formalism more flexible and prepared for more advanced applications. 
An important relation between the four types of coefficients, 

AA' BB' , An \ 
Sab = ffa ffb EABEA'B', (40) 
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can be rewritten in several equivalent ways, each of them revealing another aspect of 
their mutual relations: 

3abff a AA'S b BB' = £AB£A'B', (41) 

9 AA'fla = £A EA' , (42) 

AA' b b 

3a 9 AA' = ffa , (43) 

BA' , BA' B , AA \ 

ffaAA'ffb +3bAA'3a = E A # ab . (44) 



Eq. (|44[) is the 2-spinor form of the anticommutator known from the algebra of Dirac 
matrices. The above formulas involve boldface (numerical) indices but can be regarded 
as components of appropriate abstract-index relations, say, 

AA' BB' /A r\ 

9ab=9a gb £ ab £ a>b>, (45) 

and the like. Infeld-van der Waerden tensors g a AA , of a mixed world-vector and 2- 
spinor type, play the roles of the isomorphisms allowing to identify world- vector indices 
with pairs of the spinor ones: X a = ga AA Xaa'i an< i so forth with tensors of any rank. 

Let us note that, in general, ga AA may possess a nontrivial dependence on the 
point x of an appropriate Lorentzian manifold. The identification X a = Xaa' implicitly 
assumes that in all differential equations one encounters in a given theory one can freely 
commute ga AA with derivatives. The latter means that one works with covariant 
derivatives satisfying V a <?fc C ' C ' = (plus analogous equations obtained by raising and 
lowering the indices). In most applications to Minkowski space we do not have to worry 
about such subtleties. 

In the Penrose formalism no special role is given to Infeld-van der Waerden tensors 
— they are just regarded as one of the possible instances of the abstract-index rule 
g a b = ga BB ■ However, it is sometimes useful to work with Infeld-van der Waerden 
tensors occurring explicitly, if we need to translate some formulas into their more 
standard shapes. As an exercise, let us show several equivalent abstract-index versions 
of ([451: 

CC' DD' c d c d CC' DD' 

g a b = ga gb scd^cd' = ga gb g C d = ga gb occdd' = ga gb gcDv 

c d c d C C D D' 

~ gAA' gBB' g c d = gAA' gBB' gcC'd = £ A £A' £ B £Bi ZCD^C'D' 

CD' DC CD' DC' 

~ gAB' gBA' £CD £ C' D' = gAB' gBA' g c d = £ AB £ A'B'- 



Flexibility of switching between forms is indeed immense. We shall use this freedom in 
derivation of generators of certain important representations of SL(2,C). 



2.3 Tetrad and diad notation 

Another notational element we will need is based on the following observation. Let us 
note that, by definition of the relation between abstract and numerical indices, any 
world- vector satisfies 



X a = g a b X b = g a b X h = g a a X + g a l X x + g a 2 X 2 + g a 3 X 3 , (46) 
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so that the four world-vectors g a b play a role of a basis. Actually, this is a Minkowski 
tetrad, since 

b ac b ac be i An\ 

9a 9 = 5a 9 =9 ■ (47) 

Analogously, 

ab a cb a cb cd a b /ao\ 

9 = 9c 9 = gc 9 = g 9c 9d (48) 

a b a b a b a b / Ari \ 

= 90 90 - Si 31 - 32 32 - 33 33 (49) 

(compare Eq. (|31|) ~), A similar construction can be performed with 2-spinors, starting 
with 

<t>A = £A B <Pb = £A B 4>B, (50) 

showing that e^ B is a diad of spinorial basis vectors. In particular, since e 01 = 1 we 
find ea°£ A1 = £a°£ A1 = 1, so that o A = £a° i M = ^A 1 form a spin-frame. 



2.4 Simplest spinor representations of SL(2,C) and their generators 
SL(2,C) transformations act by <p A -> <p A = A A B (p B , cb A -» 4> A = ^/l -1 -- 4 



s 



A-I'^Vb', -> = <t> B ' A' 1 B , A ' . Here A A B and /l^ 5 ' 



the two inequivalent 2-dimensional representations of A G SL(2,C). Let us note that 
transformation properties of lower- and upper-index spinors do not have to be sepa- 
rately postulated. Indeed, for SL(2,C) transformations we have A A A B D eq d — e AB- 
Raising appropriate indices and employing antisymmetry of the epsilon, we can trans- 
form this formula as follows —A A c A B (j — e A B , which shows that —A B q = A~^q B . 

Representations (5, Jj) (in Minkowski space), (Jj,0) (unprimed spinors), and (0, ^) 
(primed spinors) of an element A £ SL(2,C) are linked to the Lie-algebra of generators 
by: A r s = ex_p^y ab a abr s , A R S = exp ^y ab a abR s , and A R , S ' = exp =£y ah a abR < S ' . 
y ab is the antisymmetric tensor whose six independent components correspond to 
boosts (t/ ™) and rotations (j/ mn ), n, m = 1,2,3. Accordingly, the generators can be 
obtained from 



dA 

°ab- = 1 



dy 



ab 



(51) 



Differentiating at y = both sides of A A Ag ecd = £ AB we obtain o ab AB = 
a abBA- Similarly, differentiating A r s = A R S A R i S , we find 

s S S' , S S / cn\ 

o a br = o- abR e R > +e R a ahR i . [bz) 

Components tr a i,r S of the Minkowski-space generators can be collected into six matri- 
ces, the ones corresponding to boosts being symmetric, as opposed to the antisymmetric 
generators of rotations. Their manifestly covariant abstract-index form, 

°abr S = i (gar9b S ~ 9br9a S ) (53) 

can be rewritten in all the possible equivalent ways by means of the tricks I have 
described above. In particular, we find 

°abr S = i {g a RR'9b SS ~ 9bRR'ga SS )• (54) 
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A A Ft R A 

Since any symmetric spinor (j>AB — 4>BA satisfies (j>A = e 4*AB — e 4>BA — 
AB 



e 4>AB = 0, comparing (J52j| with ((54j) we get 



S */ SX' SX\ ,<-r\ 

°abR = ^(9aRX'9b ~9bRX'9a ), (55) 
°abR/ S = ^(9aXR'9b XS ~ 9bXR'9a XS )• (56) 



Alternatively, after some simplifications, 



G abR S = ^£A'B'(£AR£B S +£br£A S ), (57) 
°"ahi?' 5 = ^£Ab(za>R>£b> S +£B'R'£A' S )- (58) 

We will encounter several types of spinor fields, transforming non-unitarily (unitary 
representations will appear later) and denoted as follows 

MAi...A n A' 1 ...A' m ( x ) = A At Bl ■ ■ ■^A' m Bm <l>B 1 ...B n B' 1 ...B' m ( A ^ lx ), 

(fields on Minkowski space), (59) 

MA 1 ...A n A' 1 ...A' m (p) = A A t Bl ■ ■ ■ A A' m Bm <t>B 1 ...B n B' 1 ...B' m ( A ~ 1 P^ 

(fields on mass-m hyperboloid). (60) 

A~ x p denotes the spacelike part ((A~ 1 p) 1 , (A~ 1 p) 2 , (yl _1 p) 3 ). Recall that Ap a = 
A a b p b , A P a = p bA ~ 1 b a , so that A~ 1 p a =p b A b a . 
Generators of these representations, 

J Si B^ r B\ B m i B\ B^ ( n-\ \ 

abA ± ■ ■ ■ A' m = L ab E Al ■■■£A' m +&abA 1 ■ ■ ■ A' m , (61) 

split into orbital parts L ab and the spin parts 



Bi B' . d . Bi a B' m 

°abA 1 ---AL m = l -Q^b A M A A 



y ab =Q 



(62) 



The explicit forms of the orbital parts, 

Lab<t>B 1 ...B 7l B' 1 ...B' m ( x ) = < J ab rS Xrds4)B 1 ...B rl B' l ...B' m { x ) (63) 

'dpi' 



Lab<f>B 1 ...B n B' l ...B' m {p) = ^^°ab ri Pr-^-r<f>B 1 ...B rl B' l ...B' m {p), (64) 



j = l 



in general, will not be very important if we define spins by means of Pauli-Lubanski 
vectors. 
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2.5 Bispinors and 2-spinors 

Simplest unitary representations of inhomogeneous SL(2,C) are most naturally intro- 
duced at the level of the Dirac equation, i.e. by means of bispinors. A bispinor is a 
direct sum of primed and unprimed 2-spinors. We will use the following abstract-index 

A ( V \ 

convention: ip = I A i I is a bispinor, Dirac's gamma matrices are represented by 
7,a B = ( Jq \ lq \.)= ^ ( °s g *{'), (65) 

\fqA' JqA' / \-9q A' J 

and contractions of the bispinor indices are defined by a"' A "'/3 a = a"' A "' ft A + 
....A'..., 



a P. ..A'...- Generators of the bispinor (h,0) © (0, A) representation, 

A A B A A B '\_(A A B 
A A . B A A , B ') \ A A ,# 



(66) 



satisfy 



B / O r s A B (JrsA B \ ( °~rs A B . 

^Sk = I g g, 1=1 B' 1 ' ( 67 ) 

a rsA ' Ors A i J \ U (5rs A' 



Product of two gamma matrices, 

„, B C r, ( 9qA B 9r° B' 

7gA 7''B = —2 * 



V 



ffg A'9rB 



9qr£ A ^ — 2i<TqrA 



C' C 
9qr£A' -2ia qrA ' 



c • c 

<?qr£A — 2icr gr A , 



implies 



B C B C C 

7gA 7rB + 7rA 7gB = 2g gr e A , (68) 

B C B C C 

7?a 7rB - 7»'A 7gB = -4ia qr A ■ (69) 

It is interesting that only at the abstract-index level one can realize that, as opposed 
to what one usually reads in relativistic quantum mechanics textbooks, 7oa B cannot be 
identified with the matrix /3 occurring in the standard- notation formulas ipip = ip^ /3ip 
and j a = ijjja'ip = ip ftjaip- Indeed, in the abstract-index notation the Dirac current 
roads 7 

V2g a AA ' {^ A ^A> + *Pa4a) = V2(^ A ^ A , + 4> A 4a) (70) 



Ja = 



In particular, 



f\> 5 7ob c 

, AA ' o Wv-a 
jo^'il^ 



^.toi^. if :;m- (72) 
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Replacing Infeld-van der Waerden tensors by matrices of Infeld-van der Waerden sym- 
bols we find that the matrices corresponding to /3 and 70 are indeed the same. However, 
a glimpse at the indices shows that these are matrices of maps that act in different 
linear spaces and in addition transform under SL(2,C) according to inequivalent rep- 
resentations (roughly speaking, /? A / B is a scalar, and not a timelike component of a 
world- vector) . This is an example of a relativistic formal subtlety that can be appre- 
ciated only at the abstract-index 2-spinor level. In the following sections I will discuss 
more such subtleties, clarifying certain controversies about relativistic properties of 
qubits. 

3 Qubits associated with first-quantized Dirac equation 

Although 2-spinors and qubits are vectors from C 2 , not all qubits are 2-spinors, and not 
all 2-spinors are qubits. The reason is simple: 2-spinors correspond to finite-dimensional 
non-unitary representations of SL(2,C), but qubits — by definition — have to transform 
unitarily. Qubit representations of inhomogeneous SL(2,C) are infinite dimensional. 
Physically, the infinity of dimension means that relativistic qubits are spinor fields of 
some sort. These new types of spinors are directly related to 2-spinors by a kind of 
duality. The duality will be discussed at the end of this section, but before we do that 
we have to analyze some important preliminaries. 

3.1 Relativistic generalization of spin 

One expects that relativistic qubits are related to spin. A conceptual difficulty one 
immediately encounters in this context is that there is no generally accepted definition 
of relativistic spin operator. It is not a priori evident that the discussed splitting of 
generators into 'orbital' and 'spin' parts corresponds, physically, to orbital angular 
momentum and spin. A simple illustration of the difficulty is the following. Let us take 
the Dirac Hamiltonian (units with c = 1, % = 1) 

H = a- P + mf3 (73) 

in the non-manifestly-covariant formulation. The generator of rotations reads 

J = x x P + s, (74) 

where s = 5 ^ ^ .It seems natural to identify spin with s. However, since [H, s] ^= 

the projection of s on a unit vector a does not lead to a well defined quantum number, 
unless a is an operator parallel to momentum P = — iV. The projection of spin on 
momentum P ■ s is a differential operator of first order, commuting with H. Knowledge 
of its eigenvectors is sufficient for vast majority of applications, but there are exceptions. 
Formulation of a relativistic generalization of the Bell inequality for electrons requires 
projections of spin on four different directions, so the knowledge of helicity states is 
not enough. 

A possible solution is suggested by analogous problems with velocity of a free Dirac 
particle: dx/dt = v = —i[x, H] = ct, [H, v] 7^ 0. A way out is to start with the position 
operator 



q = n + xn+ + n-xii- 



(75) 
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where 77± project on states of positive or negative energy. Then, dQ/dt = V = 
—i[Q,H] = P/H, [H,V\ = 0. Q has properties of relativistic center of mass, and 
V is the center-of-mass velocity. Applying this reasoning to spin, we arrive at 

S = n+sn+ + 77_s77_ = J Q x P, (76) 
J = xxP + s = QxP + S. (77) 

Now [S, H] = 0, and a ■ S is a natural projection of spin with no restriction on 
the direction of a. Apparently, the first application of the above idea to relativistic 
qubits (relativistic Einstein-Podolsky-Rosen-Bohm experiment) can be found in the 
unpublished preprint of mine from 1984 [3]. The resulting spin operator depends on 
momentum in a very complicated way and does not satisfy the rotation Lie algebra 
typical of angular momentum. If a = P/\P\ then a-S — as — a-J is the helicity 
operator. If a is a unit vector perpendicular to P then the eigenvalues of a ■ S tend 
to with p — > oo (in momentum space), or with m — > 0. The latter property can be 
understood as a combined effect of two classical relativistic phenomena: the Lorentz 
flattening of the particle, and the M0ller shift of the center of mass [231124] . 

Another, manifestly covariant approach is to start with first-order differential op- 
erator known as the Pauli-Lubanski vector 

Ti/ a Y r> * ji/r a k Y r> * a b Y 1 abed j-, Y /t \ 

W x = Pb M x = n <r x = 2 e p b°c<n ■ (78) 

Its projection W° = P ■ s on the 'time' direction turns out to be the helicity times 
\P\. Moreover, W a commutes with 4-momentum and W a W a is a Casimir operator of 
the Poincare group. In irreducible representations W a W a = —m 2 j(j + 1) where j is 
the dimension of the associated representation of su(2). This is why another popular 
definition of spin is w a = W a /m. Its spacelike component w equals s for p — (again, 
in momentum space). If t a is a constant world- vector then eigenvalues of t a W a tend to 
infinity with p — > oo (this is clear since in momentum space W° = p-s). It is interesting 
that S and W are proportional to each other, S = W/H, and thus possess the same 
eigenvectors. From the point of view of the Bell inequality, say, they are equivalent (cf. 
the analysis of this point given in |12II25| ). The components of W a do not satisfy the 
angular momentum Lie algebra. 

The simplest way of deriving an explicit form of W 1 -^ 1 is based on the following 
abstract-index identity: If F ab = — F ba then (cf. Eq. (3.4.21) in [7]) 



Fab = 4>AB£A'B' + ZAB^A'B', (79) 

4>ab = \f A X'b X \ (80) 

IpA'B' = \f X A' X B>, (81) 

*F a b = -i<t>AB£A'B' +ie AB^A'B'- (82) 
Comparison with (|57|) . (|58p implies 

*&abR S = ^£A'B'(£AR£B S +£BB.£A S ) (83) 

= ^(9aRX'9b SX -9bRX'9a SX ), (84) 

* a abR' S = -^£Ab(£A'R'£B' S + £B'R'£A' S ) (85) 

= -\(9aXW9b XS ' - 9bXR>9a XS ')- (86) 
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It follows that 



W aR s = \{g aR x>P sx ' - P R x>ga SX '), (87) 

WW 5 ' = -\{9aXR>P XS ' -PxR>9a XS '), 



Projections of the Pauli-Lubanski vector on various directions will be discussed in a 
separate section. 

Yet another possibility is to begin with the Newton- Wigner position operator 
QnWi which differs from Q by a term commuting with H The components of 
Qnw commute, as opposed to those of Q. One obtains Sjyw = J ~ Qnw x P whose 
components satisfy the angular momentum Lie algebra. It can be shown [2T that this 
is the only axial vector operator linear in the Pauli-Lubanski vector, satisfying su(2), 
and transforming under rotations as a 3-dimensional vector. The above properties look 
natural. Spin based on SpfW has been extensively studied in the context of relativistic 
qubits and Bell inequalities by Rembielihski and his group 28,29.301131]. 

One can argue, however, that in classical mechanics of spinning particles one does 
not arrive at commuting center-of-mass coordinates. This happens whenever one im- 
poses mass or spin constraints, which is the case here, and replaces Poisson brackets by 
Dirac brackets according to the rules of constrained dynamics. The position variable 
that commutes with respect to the unconstrained bracket becomes noncommutative 
with respect to the constrained one [321133] . and this is the bracket that should be 
employed if spin is 1/2 or mass is ?n e i ec tron- What is interesting, the resulting algebra 
is analogous to this of Q and not to Qnw- This is not so surprising if one realizes 
that position generates shifts of momentum. Mass hyperboloid in not a flat manifold 
so translations do not commute. If we relax the mass constraint, the momentum space 
becomes the Minkowski space, which is flat, so positions should commute, as it indeed 
happens in off-shell formalisms in quantum mechanics. 

Classification of all spin operators linear in momentum can be found in [34,35 , 
and different spins are used in practical applications to exact solutions of the Dirac 
equation coupled to electromagnetic fields. From a formal point of view all the possible 
spin-like observables are possible candidates for quantum yes-no observables. Which of 
them are actually measured in experiments depends on experimental procedures. 



W aR S 



3.2 Beyond helicity: Two-spinor approach to polarization operators based on the 
Pauli-Lubanski vector 

Generators of four-translations are defined by 

e*^ p 'Mx) =M*-V) = e~ yada Mx), (90) 
so P a — id a - The Pauli-Lubanski vector thus satisfies 

^ (<laHX>p SX ' -PUX'9a SX ' ). r - 

-g a XR'P XS +PXR'ga XS . 



W aR s (±p)=±W aR s (p) 

(91) 
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In all the formulas we assume that p a is future-pointing, i.e. po > 0. Now, let us 
consider any symmetric spinor cj>AB- Using 4>a A — 00° + 01 1 = 0we find 



0i° 0i X -A 



^1 



AS 



(92) 



so that the eigenvalues A' ' of <j>A B are given by 



A (±) = ±V-^xy0 xy /2. (93) 

Applying this formula to VKr S (£,p) = t a (p)W / Q R S (p), where t a (p) is an arbitrary field 
of (real) world- vectors, we conclude that the eigenvalues of such a general projection 
of the Pauli-Lubanski vector are given by 

A (±) (t,p) = ±^-W XY (t,p)WXY(t,p)/2 = ±y/-Wx>y(t, P )WX'r'(t,p)/2 

= ±±V(t-p) 2 -tV. (94) 

Here t ■ p = Pat a {p), t 2 = t a {p)t a (p), p 2 = p a p a = m . Helicity, corresponding to 
*°(p) — (1/IpIj 0, 0, 0), has eigenvalues 

\ {± \t,p) = i^pg/p 2 - m 2 /p 2 = ±|. (95) 

If t a (p) = p a then A'^^p) = 0, since P a W a = in any representation. The latter 
incidentally shows that projections of the Pauli-Lubanski vector possess a gauge free- 
dom: t a W a = (t a + 6P a )W a , for any 9. The same concerns the eigenvalues: \ t ' ± \t,p) = 
\( ±S) (t + 0p,p). 

For a particle at rest, p a = (m,0), the Pauli-Lubanski vector reduces to non- 
relativistic spin multiplied by mass: W a = m(0, s). The non-relativistic observable 
as— ■ W can be covariantly written as t a W a , where 

t a = (t°, a/m) = (0, a/m) + (t° , 0), (96) 



the last term being proportional to p a . This is the simplest explanation of the gauge 
freedom inherently present in the Pauli-Lubanski vector. The fact that it is only the 
spacelike part of t a that counts in the definition of t a W a does not mean that t a itself 
should be spacelike. A spacelike t a when shifted by a timelike 8p a can become timelike 
or null (or remain spacelike). 

Of particular interest turns out to be the null case: t 2 = 0, t a p a = 1. Indeed, first of 
all the corresponding eigenvalues are identical to those of the helicity, A^'(t,p) = ±4. 
Let t a = LJ a (p) = uj A (p)w A ' '(p), and define 

n A (p) = p AB w B ,(p). (97) 



The trace-reversal formula (Eq. (3.4.13) in [7]) 



2 2 
AB' BA' AA' BB' m AB A' B' a b m ab /na \ 

p p = p p - — e e =p p — g , (98) 
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implies 

A, \-A' i \ AB' BA' i \- , \ 

TT (p)7T (P)=P P LU B {p)u) B '{p) 

m 2 i 

= p a p b uj B (p)ui B '(p) tt^ (P)u (P). (99) 

p b t b =l 

and we arrive at the following important decomposition of p a into a combination of 
two null directions: 

' m 2 ' 

p a = tt A (p)tt A (p) + —uj A (p)ui A (p). (100) 

Formula (|100|) implies ([97]) • Indeed, 

2 

AB' - ! \ / A/ n-B'/ n . m A, \-B' , A - / \ A / \ 
P LU B ,(p)=\tT (p)tT (P) + — W (p)w (p) 1 W B , (p) = 7T (p). 

It follows that uja (p) supplemented by (| 100p determines its spin-partner ir A (p) uniquely. 
Moreover, since 

' m 2 

PAC'K A (p) = PAC'P AB u B >(p) = — u C '(p), (101) 

the knowledge of iv A (p) and (|100p determines kM(p) uniquely, unless m = 0. 

As opposed to the special case m = 0, where n A (p) is defined by its flagpole p a up 
to a phase (U(l) internal symmetry), the internal symmetry group is here U(2). The 
pair (tt A (p), -^uj A (p)) can be replaced by (ir A (p), -^=a> A (p)), where 

\^5 U {p) ) Vt(p) 5 (pV {p) ) 

and the matrix is unitary. Hughston showed 36 that the existence of such splitting-of- 
p a ambiguities may be the actual geometric reason for internal symmetries occurring in 
gauge theories. From our point of view it is more important that a similar mechanism 
is responsible for the structure of unitary representations of the Poincare group |22| . 



3.3 The simplest case: Projection on the null direction t a = uj a (p) 

Since uj a p a = lu A uj A tta^a' = 1 independently of the value of m, we can assume — 
without loss of generality — spin-frame normalization oja 7vA = 1- Projections on the 
null direction, 

W(uj,p) A B = i (n A w B + ^tt S ) , (103) 
W(u,p) A , B ' = -\(^a^ B ' +u A 'K B '), (104) 
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look the same for both m = and m 7^ 0. Also spin eigenvalue problems become in 
the null formalism particularly simple: 

B 1 

W(uj,p) A oj b = -u A , (105) 

W(oj,p) A , B 'n B , = (106) 

W^,p) A S 7r B = -~7T A , (107) 

Wfw.pU^W = -iffl A /. (108) 

Spectral projectors associated with X^\ui,p) = ±4, 

n {+ \oj,p) A B =uj a it B , (109) 
I^W)^ = -7T^ B (110) 

n^\u,p) A , B ' =-* a ,q b ', (in) 

TT^.p)^' = cW\ (112) 
trivially satisfy 2-spinor resolutions of unity (|29|l , 

ZZ" (+) (w,p)^ fl + 7T ( ~ ) (w,p)a' B = ^att B - tta^ B = (113) 
il(-)(w,p) A / B, +irW(a;,p) A /- B ' =^7f B '-7f A ^ B ' =e A , B '. (114) 

At the bispinor level the projectors read 

n(+)< \ ^ ( n<- + \u,p)A B \ /w A 7r B \ 

77 (WiP)a =1 ^w^'M J- :m; 

r(-)( r , n \* - f ^ (_) (w,p)A B \ = f-v A u B 

n(-\u, P ) A , B ' J \ qa,* 1 



3.4 Wave functions associated with W(oj,p) B 

In order to introduce wave functions corresponding to W(u),p)k B we first have to ex- 
pand solutions of Dirac equation in terms of appropriate eigen-bispinors. A systematic 
procedure of deriving the bispinors can be based on 'spin-energy' projectors. 
The free Dirac equation 

£>aVb = 0, D B = id r -yrk B - me B (117) 

is implicitly an orthogonality condition for 'sign-of-energy' projectors II±k B (p) con- 
structed as follows. One begins with 

D B e TiP-x = £> ±A B(p) e =FV-* ; ( 118) 

where 

pi b / \ , r b „, B ( -me A B , ±V2p A B \ . 
D±k (P) = ±p Irk - msk = r- g gi I ■ (119) 

V=FV2p A' , -raE A , a ) 
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The orthogonality conditions 

D±, B ( P )D TB c (p) = 0, (120) 

D± k B (p)D ±B c (p) = -2mD ±A c (p), (121) 

imply that the possible eigenvalues of -D± A B (p) are and — 2m. Accordingly, the pro- 
jectors we are interested in read 

77 ±A B (p) = D ±A B (p)/(-2m) = \ [ *f T ^ P f ) . (122) 

2 V ± ™ p A ' £A ' J 

Multiplying (|122[) by (j 1 15f> or 1)11611 we obtain 'spin-energy' projectors, 

4+ Wp) = 27 ± aW +) . Vp) = \ ( c T ^ A f ) , (123) 

n y ± \ {oj,p) = n ±k (p)m >b (w,p) = - m _ c _ C ' • ( 124 ) 



o I i rra - o — — (. 



The completeness relation is 



2 i7i S V(^ P )= £ A B =( e f J^). (125) 

It is interesting and important for later applications that (|123p and ()124p allow us 
to define basis bispinors that can be used also in the massless case. The trick is the 

following. Let <j> A (p) = 2 ( " F\ ] ■ The eigen-bispinors 
V n A' (p) J 

^\(u, P ) = i7£ + > A Vp)<M P ) = ( ±J r^ f A(p) ) , (126) 

V n A> (p) / 

4> { ± \(u,p) = n ( ±\\u,p)Mp) = ( T ^o^(p) ) . ( 127 ) 

are well defined also for m = 0, and satisfy 

W(w,p)i , V£ ±) b(">p) = ±^</>i ± \(u,p), (128) 

Bi^h^p)^*"* = -2m4>¥\(u>,p)e* ip - x , (129) 

D A B 4 s ') B ( W ,p)e* = 0. (130) 

A general solution of free Dirac equation can be finally written as 



il>i(x) = / dp{i>-k{p)e tp ' x +ip+k(p)e +zp ' x ^J 

J dp(<^ (s) A^,p)/( S ,p)e-^ a; + ^ ) A (a;,p)^F^) e + ^ a: ). (131) 



s=± ■ 
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Here dp = d z p/[{2-K) z 2y/p^ : + m 2 ] is the invariant measure on mass-m hyperboloid, 
and complex conjugation and the minus sign in g(—s,p) are convenient for later appli- 
cations. By definition, 

mp) = ( ^ Ap) ) h + ,p) + (_ :^LV(-, P )- d32) 



^«=("?xr«^H#-W^ (i33) 

Wave functions f{s,p) and g(s,p) are simultaneously scalar fields and spinors of a new 
type, as we will see shortly. They can be extracted from ip±k(p) in a simple way: 

f(+,p) =Q A >(p)4>+ A '(p), (134) 

/(-,p) = o, A (p)i/, + A (p), (135) 
A' 



ff(-,p) = <*>A>(p)il>- (p), (136) 
= wa(p)V>- A (p)- (137) 

For future reference let us explicitly write the two components of the m = case: 



1> A (x) = J dfc7r j4 (fe)(/(-,fe)e-* x +g (+,k)e ik - x y (138) 
iM*) = f dfc^(fe)(/(+,fe)e- lfc - a; +^rfc)e <fc - x ). (139) 



Massless 4-momenta are denoted here by A; a and dA; = d 3 fc/[(27r) 3 2|fe|] is the invariant 
measure on the light cone. 

This is the right moment to explain the issue of unitary representations of the 
(covering space of the) Poincare group. I will not follow the usual Wigner-Mackey 
procedure [1811 19] of induction from little groups. Instead, I will show by means of 
2-spinor techniques that unitary representations of the group are encoded in shapes 
of solutions of relativistic wave equations. The latter statement is in itself not very 
original since links between unitary representations and relativistic wave equations 
were discussed already in [20] (for a modern analysis cf. [2T]). Nevertheless, the 2-spinor 
tricks I will use do not seem to be widely known and, apparently, were introduced for 
the first time in 22 . 



3.5 Duality between active and passive SL(2,C) transformations: cj-spinors 
Let us now consider the transformed solution 

= J2 I dpV(^ ) B^,p)/( S ,p) e - ip - A " x + ^ ) B(^,p)^F^) e +lp ' A " a; ) 

s J 

Employing p- A~ x = Ap-x, changing variables in integrals, keeping in mind that dp = 
d(Ap), and finally expressing the transformed solution again in the basis <fi±\(u),p), 
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we get 



AM*) = Y1 f dp(A k *<t> ( °\{u 1 A- 1 p)f{ S ,A- 1 p)e-^ x 

s •* 

+ V<^W, A- x p)g{-s,A-lp)e +i P *) 
= E/ dp{^\{w,p)Af{s,p)e-^ x + ^\(u;,p)Ag(~s,p)e + ^ 



Comparing appropriate terms we find 



Att A '( P ) 



7f A '( P ) 



f(+,A- 1 p) + 
Af(+,p) + 



An A (p) 
n A (p) \ 



V2 



co A (p) 



(140) 

4f(-,p), (141) 



where An A (p) = A A B n g(A~ 1 p), Auj A (p) — A A B ujg(A~ 1 p), and analogously with 
the primed spin-lrames and complex-conjugated anti-particle wave functions Ag(—s,p). 
Denoting 



we obtain 



\tp+l {p), 



M±a(p) = C/(ApU 8 ^±e(^" 1 p), 




(142) 



(143) 



( H±o(p) \ _ ( u A ,(p)An A '{p) -^u A ,(p)Au A '(p)\ / 1 / )±0 (A- l p)\ 
\M±1{P)J ~ \^u, A (p)Au: A (jp) u; A (p)An A (p) J {^(A^p) J ' 

(144) 

Notice that new a type of binary indices has been introduced: A, B = They 
correspond to local SU(2) spinors. (Local since spinor fields taken at different points 
of the mass-m hyperboloid transform by different SU(2) transformations.) 

It is obvious from the construction that the latter passive [7] transformation of wave 
functions is equivalent to the bispinor-field acive transformation AiPa(x) = A^tp^{A~ 1 x) 
of solutions of the Dirac equation. One can explicitly check 22 that the matrix in (1144 [ 
is unitary and has unit determinant, and that the map f(s,p) — > Af(s,p) is a rep- 
resentation of SL(2,C). The choice of f(s,p) is completely arbitrary and unrelated to 
the Dirac equation itself, but the equation is encoded in the form of the matrix occur- 
ring in (|144p (via its implicit dependence on 'sign-of-energy' projectors). Translations 
x a — > x a + y a are represented at the level of wave functions by f(s,p) — > e zp ' y f(s,p). 
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The map f(s,p) — > e lpy Af(s,p) is thus nothing else but the unitary spin-1/2, mass-m 
representation of inhomogeneous SL(2,C). The massless limit m = is trivially ob- 
tained from (|144[) and shows that the representation splits for m = into direct sum 
of two representations. 

The above properties are well known in the context of induced representations of 
inhomogeneous SL(2,C) [18] : The case m > corresponds to the little group SU(2) 
of p a — (m,0); for m — the little group is E(2) and its irreducible discrete-spin 
representations are 1-dimensional. The formulation I have presented in this section 
depends on two crucial technical elements, the decomposition (|100|) and the choice of 
t a = ui a , but does not use the idea of induction from little groups. 

Local SU(2) spinors associated with projections of W a on t a are termed the t- 
spinors. A link between general t-spinors and the oj-spinors is given by an SU(2) trans- 
formation whose explicit form can be found in |22| . For a non-null t a the transformation 
is much more cumbersome than the case of null uj a , but choosing t a = (l/|p|,0) we 
reconstruct the standard helicity formalism [37] . The direct sum (|142|l of particle and 
anti-particle representations can be called an o;-bispinor. 



3.6 Pauli-Lubanski vector associated with (|144p 



It pays to compute explicitly the Pauli-Lubanski vector corresponding to (|144l) . 



, n B b* ■ d 
Wa{P)A = P l 



u, A ,(p)AK A '(p) -^a A ,{p)A^ (p) 



The projection 



dy<* yj=u; A (p)Au; A (p) u A (p)An A (p) 
' Cj R '(p)Wa(p)R> S 'n s ,(p) -^ R ' (p)W a (p)w S 'oJ s ,(p)\ 

^ R {p)W a {p) R S U S {p) CO R (p)W a (p) R S 7V S (p) J 

7r^O)7fA'(p) - t w a(p)^A'(p) -mV2K A (p)u A >(p) ) 

-mV2uj A (p)n A , (p) ^-uj A (p)ui A , (p) - tt A {p)k a , (p) J 

(145) 



u A (p)u A '(p)W a (p) A B = W{u,p) A B 
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is the same for both m — and m > 0. For one can employ 

Al \ - A' I N TJ/ / N I 



U A (p)n A ' \p)Wa{p) A B 



to define relativistically invariant analogs of the remaining two Pauli matrices, 
The massless case, 



W a (k) A L 



1 

-1 



ha 



(146) 

(147) 
(148) 

(149) 



explains why projections of the Pauli-Lubanski vectors cannot violate the Bell inequal- 
ity if m = 0: All such observables commute with one another. The ultrarelativistic 
limit of massive particles is more complicated 12,28 . 
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3.7 Peres-Scudo-Terno problem and an w-spinor way to circumvent it 

True relativistic qubits are, by definition, fields that transform unitarily under inho- 
mogeneous SL(2,C), i.e. the t-spinors, or the i-bispinors. Let us concentrate on the 
t-spinor case. The corresponding Hilbert space is equipped with the scalar product 



The expression p(s,p; s ,p') = f(s,p)f(s',p') is an example of a pure-state one-particle 
density matrix of a relativistic qubit. The qubit is characterized by two quantum num- 
bers — spin s and momentum p. Peres, Scudo and Terno posed the following problem 
[5]: Assume that our detectors do not distinguish between different p, but detect the 
sign s in /(s, p). Can we use this s as a quantum bit? The probability of finding a given 
s is given by J dp p(s,p; s, p). It is therefore justified to introduce the reduced spin den- 
sity matrix p(s,s) = J dp p(s , p; s' , p) . The problem is that the entropy of p(s,s') is 
not relativistically invariant. Indeed, let f(s,p) = F(s)G(p), where J dp\G(p)\ 2 = 1, 
\F(s)\ 2 = 1, and p(s, a') = F(s)F(s'). The wave function F(s) is a pure state so the 
entropy of p(s, s') is zero. Now, the matrix U (A, depends on momenta in a highly 

nontrivial way, and its form is even more complicated if one works with helicity or other 
t-spinors. In consequence, Af(s,p) does not separate into a product AF(s)AG(p), and 
thus involves a nontrivial entanglement between s and p. But then it is well known 
that the reduced density matrix Ap(s, s') is mixed: Its entropy is nonzero. The result 
has attracted attention of many authors and various variations on the theme can be 
found in the literature. With the exception of 9,38 apparently all the authors work 
in the helicity formalism. 

But we know that qubits can be defined in terms of all t-spinors — there is nothing 
fundamental about helicity, a t-spinor associated with direction parallel to some t a — 
(1,0,0,0). In particular, let us take an arbitrary null direction t a = t A t A . If p 2 = 
m 2 > then p ■ t is nonzero for all p a . The null vector ui a = t A t A /(p ■ t) satisfies 
p-ui = 1, so oj A {p) — t A /yjp ■ t can be used to split 4-momentump° according to (|100p . 
and the whole construction of w-spinors can be repeated. Now, let A^ B be an arbitrary 
SL(2,C) transformation. It is known (Sec. 3.6 in [7]) that any SL(2,C) transformation 
possess at least one eigenvector, A^ B rg = Ar^, A = \\\e lv , so let us take this in 
our definition of t a . Then, 




s=± 



(150) 



Aa B u b {A 1 p) 



e- lip u A {p), 



(151) 



A A B 7v B (A- 1 p) 



VWVg a a b a-i b c pc b 'a b , c 't c , 




(152) 



The transformation matrix becomes 



( 



u A '(p)An A '(p) -^A'(p)AoJ A '(p) 
^ A (p)Au 1 A {p) u> A (p)An A (p) 



) 
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_/ e^u> A ,(p)^ A '(p) - e -^^Cj AI (p)cJ A \p)\ = f e w o \ 

\e^^ A (p)u A {p) e-^u A (p)n A (p) ) {0e-^)' [b6) 

which is independent of p. Actually, the transformation matrix is simply an ordinary 
rotation by 2ip which, of course, cannot change any entanglement. The reduced density 
matrix of Peres, Scudo and Terno has not changed its entropy, even though the SL(2,C) 
transformation was arbitrary. 

The construction I have just presented was introduced in [9], and can be regarded 
as a relativistic error correction algorithm. What it means is that it is always possible to 
adjust t a in t a W a in a way that ensures that qubits defined by this concrete observable 
do not change their Peres-Scudo- Terno entropy. This trick can be generalized to any 
number of particles that propagate in arbitrary ways, not necessarily inertially, and 
described by any states — entangled or not. Another consequence is that entanglement 
between Peres-Scudo- Terno qubits will not change, no matter what kind of motion is 
considered, if one defines qubits by means of t a which are principal null directions 
[S] of Lorentz transformations that define the motion of observers. 



4 Transformation properties of fields of spin frames for massless particles 

Let us consider a massless particle whose 4-momentum k a is a future-pointing null 
but nonzero world-vector, k a k a = 0, k a 7^ 0. We know that such a k a determines, 
up to a phase, a 2-spinor n A satisfying k a — n A it A . The latter equation can be 
understood in several ways. The approach I prefer takes the future light cone as a 
given 3-dimensional manifold on which a spinor field n A — Tv A (k) is defined, but 
treats this spinor field as a fundamental physical object related to 4-momentum by 
k a = n A {k)n A '{k). A consequence is apparently paradoxical. The spinor field trans- 
forms under SL(2,C) transformations as n A (ft) h-> An A (k) = A A B Tr B (A- 1 k). But 
then A-k A {k) A-k A i {k) — A a b A~ 1 i } c kc = k a - Accordingly, the 4-momentum does not 
change. This type of approach is consistent with the idea that the fundamental fields 
are those of 2-spinors and not the 4-vectors. An important consequence is that both 
An A (k) and n A (k) possess the same flagpole and, thus, differ at most by a phase fac- 
tor: A"K A (ft) = e~ l0 ( A < k K A (k). For any spin- frame partner 1^4 (ft) of Tt A (k), we find 
U) A (k)Air A (k) = e ~ i0{ - A ' h ^ and Q A > {k)~AK A> (ft) = e ie(A ' fe) . The latter two expressions 
have been encountered before in the diagonal elements of the matrix in H144[) . 

Now let us consider the 'orbital' part of the 6-angular momentum antisymmetric 
tensor of the massless particle: M a ^ — x a k^ — x^k a - Any antisymmetric (real) second- 
rank tensor can be written as (Eq. (3.4.20) in [7]) 

M ab = MASEA'B' + £ABP-A'B', (154) 

1 a' b' 

VAB = 2 M AA>BB>£ = HBA, (155) 

Pa'b' = ¥ab, (156) 

the relation between M a f, and h A b being expressible also by means of the generator 
oiA A B : 

Mab S T.fab' 1 / SX' SX'\ ■ ji^-ab SX' 

°abR =M -{g a RX'9b -9bRX'9a )=lM g a RX'9b 

= iM RX , sx ' =2i m S . (157) 
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Let us apply this procedure to M a {,: 

HAB = ^(^afefc - X b k a )e A B 

1 / - — \ A'B' 

1/ -A' -B' \ 

= 7}(u A TT B +U) B TT A ). (158) 

The orbital part thus turns out to be determined by two spinors: n A whose flagpole is 
the 4-momentum, and u) A — u) A (x,ir) = x AA in A . 

The bispinor (uj A ,n A i), with the characteristic bilinear dependence of uj a on x 
and 7f, is known as a twistor [Sj. Treating n A as a spinor field 7r^(fe) we obtain 
oj^ = uj A (x,k) = x j 4 J 4'7f' 4 (fe) which is a field defined on the Cartesian product of 
the Minkowski space with the momentum-space light cone. The field transforms under 
SL(2,C) transformations as follows 

u A (x,k) — > A A B uj b (A~ 1 x, A 1 k) 

= A A B A- x B c A- 1 B , c 'x C c^ B '(A- 1 k) 

-l C 



= e A C x cc ^ B \A- 1 k)A- 1 B , 

= x ac ,Att C (fe) 

= x AA >e 0{A ' k) TT A (k) 

= e i0(A ' h) u A (x,k). (159) 

The contraction 

^(a;, k)ir A (k) — x AA > ir A (k)it A (fe) = x a fc a = x ■ k (160) 
shows that for x a fe a 7^ we can define a spin- frame field 

M(fc) = Jr A (fc), (161) 
0a(X; fe) = fOA&Q. = -AA'n A '(k) ^ (162) 

satisfying simultaneously 

o A (x,k) A (k) = 1, (163) 

L A (k)T A '(k) = k a , (164) 

A A B L B (A^k) = e^'tAffc), (165) 

A A B o B (A~ x x, A~ x k) = e <e(A ' fc) o,i(ir, fe). (166) 



It is an instructive exercise to check what would have been the consequences of assuming 
that x AA itx A (k)/(x • k) is a spinor field defined on the light cone only, and not on the 
'phase space' (x,k). Then we would have to assume the transformation rule 

Ao A (k) = A A B o B {A~ X k) = e l0{A ' k) Ax AA ,n A ' \k)/{Ax ■ k). (167) 



25 



The right side of (|167|l is not proportional to o A (k) but also contains a component 
parallel to i A (k). In effect, we find 

Ao A {k) = A A B o B (A- 1 k) = e i0{A ' k) o A {k) (168) 

where, in spite of the fact that o A (k)b A (k) — 1, o A (k) — o A (k) + (j)(A, k)i A (k). The ap- 
pearance of the function cf>(A, k) = o A (k)d A (k), which is different for different explicit 
forms of the spin frame, is responsible, as we shall see later, for a gauge transforma- 
tion of the electromagnetic four-potential and has consequences for field quantization. 
This change of gauge seems to be a generic property of spin frames defined on mass-m 
hyperboloids. 

However, we will see that there are reasons to believe that Eq. (|166[) is of funda- 
mental importance for relativistic transformations of electromagnetic qubits, but one 
should not identify the point x a with an event in Minkowski space. A more natural 
interpretation arises if one treats the twistor as defined with respect to some internal 
coordinate, say R a , which is timelike (and thus R ■ k is nowhere vanishing if k a 7^ 0). 
At the level of quantized electromagnetic field the approach will lead to an intriguing 
splitting of the field in momentum space into three massless fields: one field trans- 
forming according to the unitary spin-1 representation of the Poincare group, and two 
scalars. Photon degrees of freedom will be defined with respect to the spin-1 part. 

5 Qubits associated with electromagnetic field 

Qubits associated with electromagnetic field can be discussed at both first- and second- 
quantized levels. In this section we concentrate on the first-quantized description, which 
is important since transformation properties of polarization degrees of freedom will 
remain unchanged after field quantization. 

Electromagnetic field tensor is antisymmetric, so can be decomposed in the usual 
way: F ab (x) = F AB (x)e A , B , + £ AB F a , b , (x), F a , b ,(x) = F AB {x), F AB (x) = F BA (x). 
Source-free Maxwell equations, d a F a f ) {x) = d a * F ab {x) — 0, are equivalent to the spinor 
equation 

d AA 'F AB (x) = Q = d AA 'F A , B ,(x). (169) 

In Fourier space 

F AB {x) = Jdk (F_ AB {k)e- lk x +F +AB (ky k - x y F ±AB (k)n A (k)7t A '(k) = 0. 

(170) 

Expanding F± AB (k) in terms of spin frames we find that the only symmetric solution 
of (|170|l is F± AB (k) ~ n A (k)n B (k). The field spinor can be thus written as 

F AB (x) = -J dkiv A (k)n B (k) (a(- fc)^"* + a(+, k)e lkx ^ , (171) 

F A 'B' (as) = - J dk n A , (k)-K B , (fc) (a(+, k)e- lk ' x + a(-,fe)e* x ) . (172) 

The signs and conjugations are a matter of convention. Similarity of (|171|) - (|172|) to the 
zero-mass solution of Dirac's equation (|138|) - (|139|) is not accidental and is typical of 
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massless Bargmann-Wigner equations of any spin [221411 . SL(2,C) transformations of 
(|171I) - (I172[1 can be deduced from analogous formulas found for the Dirac field, 

A A c A B D F CD (A- 1 x) 

= - J dkit A (k)7r B (k)e- 2ie(A - k) (a(- A 'A-:- ' ? " + a(+, A-^-k)e lkx ^j , 

(173) 

A A , C ' A B , D ' F C , D ,{A- X x) 

= - J dk tt a , (k)7v B , (k)e 2z0{A ' k) (a(+, A' 1 ^-^ + aC-.ii-ifcJe**'*) . 

(174) 



5.1 Four-potential A a (x) based on spin frames defined on the light cone 

Four-potential is defined by 

F ab (x) = d a A b (x) - d b A a (x). (175) 

We know that given -K A (k) we cannot find a unique spin-partner oj A (k), since for any 
function cj>(k) the new field 

Cj A {k) =cj j4 (fc) + 0(fc) 7 r A (fe) (176) 

also satisfies uj A (k)n A {k) — 1. This leads to the equivalence relation: Cj A (k) ~ ui A {k) 
iff Cj A (k) — ui A (k) is proportional to iv A (k). 

The 4-potential in a Lorenz gaugtQ, can be taken in the form (cf. [39U40] ') 

A a {x) =i J dke~ tk ' x (u A (k)n A ,(k)a(+,k) + n A (k)uj A ,{k)a(-,k)^ + c.c.(177) 

Now, if we replace uj A (k) by Cj A (k) satisfying (1176[1 . then 

) - da${x) (178) 

where 

${x) = J dk 4>(k) (a(+, k)er ' lk ' x + a(-,k)e ik ' x ^ + c.c. (179) 

is a solution of d a d a &(x) = 0. The equivalence class of spin-frames corresponds to an 
equivalence class of Lorenz-gauge potentials. 

(|173p ~ (|174p show the duality between active spinor transformations and passive 
zero-mass, spin-1 unitary transformations of wave functions. An analogous construction 
does not exactly work for A a (x). Let us have a closer look at this important subtlety 
since it has nontrivial implications for field quantization. 



Note: 'Lorenz gauge' (L. Lorenz) and 'Lorentz group' (H. A. Lorentz). 
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The problem is with changes of gauge of the form (|178|) and (|168[) : 

A a b A b (A- X x) = i J dke~' k - x e 2ie{A ' k) Cj A (k)7T A ,(k)a(+,A- 1 k) + ... 

= i J dke' lk - x e 2l0{A ' k) u A {k)n A ,{k)a(+, A _1 fe) + d a $(x) + . . . 

(180) 

In consequence, if wave functions transform according to massless spin-1 unitary rep- 
resentations of inhomogeneous SL(2,C), then 

A a (x) -> A a b A b (A~ x x) - d a ${x), (181) 

meaning that A a (x) is not a four- vector field. And the other way around, if we assume 
that the four-potential is a four-vector field, its corresponding wave functions do not 
carry massless spin-1 unitary representations. 



5.2 Null vs. Minkowski tetrads, or circular vs. linear polarizations 

The polarization world vectors occurring in (I177|l can be regarded as elements of a null 
tetrad in Minkowski space (we skip the arguments k): 
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and dually 
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(183) 



The associated Minkowski tetrad can be expressed by means of Infeld-van der Waerden 
symbols, g a BL = g a BB 5°bb', 
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and dually, g a & = g & BB'9a 
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The meaning of x a and y a can be deduced from 

x a {k)a\(k) + ?/a(fc) a 2(fc) = TO a (fe)a(+, k) + m a (k)a(—, k), (185) 
a(±,k) = -j=(ax(k)±ia 2 (k)j. (186) 

Elements m a {k), fh a (k) of the null tetrad correspond to two circular polarizations 
(since after SL(2,C) transformations they get multiplied by appropriate Wigner phase 
factors). The two elements of the Minkowski tetrad, x a (k), y a (k), define linear polar- 
izations, 

A a (x) = i J dke- ihx (g a 1 (k)ai(k) + 3a 2 (fc)a 2 (fc)) + c.c. (187) 



5.3 A a (x) based on twistor-like spin-frames 

Let R 2 = 1, Rq > 0, dR = d 3 7?/[(27r) 3 2Vl + -R 2 ], and consider spin-frames analogous 
to (fT6T > -(fT56 |> : 

LJ A {R,k)TT A (k) = 1, (188) 

7v A {k)n A ,{k) = k a , (189) 
A A B 7v B {A- 1 k) = e- t0{A - k) ir A {k), (190) 
/U Vb(^ _1 -R, ^ _1 fe) = e j6l(yl ' fe) a; A (ii, fe). (191) 

The world-vector R a has properties of four- velocity. In this sense it is analogous to the 
four-velocity parameter occurring in theories with preferred reference frame or based 
on non-standard clock synchronization schemes 31 , 42 , 43 . Another interesting analogy 
is with fields defined by Kaiser on phase-space |44[|45| — here R a would be an analogue 
of the temper vector. The four-potential 

A a (x) — i I dk dRe~ zk x 



x(^j A (R,k)-K AI (k)a(+,R,k) +-K A (k)u} A ,(R,k)a(~,R,k)j + c.c. 
= i J dkdRe~ lkx (g a 1 {R,k)a 1 (R,k) + g a 2 {R,k)a 2 {R,k)j + c.c. (192) 

transforms as follows 

A a b A b (A~ 1 x) = i J dkdRe^ lk ' A ' lx A a b ujB{R,k)n B ,{k)a(+,R,k) + ... 

= i J dkdRe^ lk ' x A A B uj B (R,A~ 1 k)A^ BI (k)a{+,R,A~ 1 k) + ... 

= i J dkdRe~ lk - x e 2i0(Ak) uj A (R,k)n A ,(k)a(+,A- 1 R,A- 1 k) + ... 

(193) 
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Note that linear polarization vectors are elements of the Minkowski tetrad transforming 
as follows: 



/ At a {R,k) 
Ax a {R,k) 
Ay a (R,k) 

\Az a (R,k), 




/ Auj A (R,k)A^ A ' (R,k)\ 
Auj A (R,k)A^ A ' (k) 
ATr A (k)A7J A ' (R, k) 
V A-K A (k)A^ A '(k) J 

t a (R, k) 

cos 2G(A, k)x a (R, k) + sin2(9(/l, k)y a (R, k) 
- sin 20{A, k)x a (R, k) + cos 2G(A, k)y a (R, k) 
z a (R, k) 



(194) 



(195) 



t a (R,k) and z a (R,k) do not appear in (|192[) . but the space-like plane spanned by 
x a (R,k) and y a (R,k) is invariant under SL(2,C). So the twistor-like spin-frames 
lead to Minkowski tetrads that span two orthogonal subspaces, both invariant un- 
der SL(2,C). This does not hold for tetrads corresponding to spin-frames that depend 
only on k. 

The passive unitary representation of the amplitudes reads 



Aa(±, R,k) = e 



±2i0{A,k) 



a(±,A~ 1 R,A~ 1 k). 



(196) 



There is no change of gauge. (Change of gauge means that the plane spanned by m a 
and ifi a is not invariant, since a contribution parallel to k a occurs. In quantum-field- 
theory terms this would mean that SL(2,C) transformations mix transverse photons 
with 'longitudinal or timelike' spin-0 'photons'; see the discussion in the next Section.) 
The field spinors 



Fab(x) 



dk dR ir a 



(k)ir B (k)(a(- 



, R 1 k)e lk ' x + a(+, R, k)e 



ik-x 



Fa'B'(x) = - / dkdRn A ,(k)n B ,(k)(a{+,R,k) e - lk ' x + a(- , R,k)e lk ' x 



(197) 



,(198) 



are identical to (|171[) - H172[I if J dRa(±, R, k) = a(±, k). Moreover, relativistic invari- 
ance of dR implies that 

Aa(±,k) = e ±2l0{A < k) J dR a (± y A- 1 R,A~ 1 k) = e ±2t6{A ' k) A^k)^) 



Let J dR aj (R, k) = aj (k) , j = 1, 2. SL(2,C) simply rotates the amplitude 2-dimensional 
vector (ai(k) , ct2(k)) by spin-1 Wigner angle 20(A, k). 

The world-vector R a plays a role of an internal degree of freedom, invisible at the 
level of field tensors. Its role can be fully appreciated only at the level of quantized 
fields |10j . Let us stress that the choice of a time-like R a is perhaps not necessary — in 
principle R a could be space-like or null, but only for a time-like R a the product R ■ k 
is, for k a ^ 0, non-vanishing, which allows us to define spin-frames globally for all k a 
and R a . The globality condition is convenient but not necessary. 
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6 Qubits associated with quantized electromagnetic field 

An automatic identification of the Pauli-Lubanski vector with polarization of a quan- 
tum field leads to a problem: W a is quadratic in elements of Poincare algebra, W a = 
Pb* M ab , and thus — in terms of number operators — W a is quadratic as well. However, 
spin observables, similarly to 4-momentum and 6-angular momentum, are expected to 
be linear in numbers of particles. The subtlety is not so visible in first quantization 
since observables are there defined at 1-particle levels. A possible way out is to treat 
the harmonic oscillator algebra as the fundamental algebraic level, and treat inhomo- 
geneous SL(2,C) as a structure derived from it. The choice of spin-fames — possessing 
appropriate transformation properties — becomes then essential for a correct relativis- 
tic formulation of the formalism. 



6.1 Manifestly covariant quantization of A a (x) based on twistor-like spin-frames 

One typically quantizes A a (x) in a way that does not guarantee the four- vector behavior 
of the field. An approach to relativistic EPR correlations of photons based on such a 
non-manifestly covariant formalism was recently discussed in |11II46| . In what follows 
I will outline an alternative but not widely known formulation based on tetrads built 
from twistor-like spin-frames [10] , 
The potential is defined by 

A a (x) = i I dk dRe~ lk x 



: (g a l (R, fe)oi(H, fc) + g a 2 {R, k)a 2 (R, k) 

+ g a 3 (R, k)a 3 (R, k) + g a °(R, k)a Q (R, fc) 1 ) + c.c. (200) 



The field £r a a (-R, fc), by definition of a Minkowski tetrad, satisfies g a & {R, k)g ah (R, fc) = 
<? ab , g aBL (R, fc)<7f, a (R, fc) = g a b- The amplitude operators are assumed to satisfy the 
ordinary harmonic-oscillator Lie algebra, 

[a a (i?,fc),a b (-R' ! fc') t ] = S ah S(R,k,R',k')I(R,k), (201) 
[a a (R, k),n h (R', fc')] = S ah 5{R, fc, R 1 , k')a a (R, fc), (202) 

[aa(R,fc) f ,n b (i?',fc')] = -5 ab 5(R,k,R',k') aa .(R,k)\ (203) 
[a a (R, fc), I(R', fc')] = [a*(R, fc) f , /(#', fc')] = [n a (R, k)\l{R', fc')] = 0, (204) 
[oa(H, k),a h (R', fc')] = [a*(R, k)\a b (R', fc') f ] = 0. (205) 

Here <5 ab , a, b — 0, 1, 2, 3, is the Kronecker delta; S(R, fc, R' , fc') = 8i(R, R )So(k, fc'), 
where the latter two distributions are SL(2,C) invariant Dirac deltas on R 2 = 1 and 
k 2 = hyperboloids, respectively. Central element I(R, fc) is a multiple of identity 
in irreducible representations of the algebra. However, I believe that a class of re- 
ducible representations discussed below is more physical 40,10,47] so I keep I(R,k) 
at this stage unspecified. How to choose the representation of the algebra becomes 
particularly important if one considers the so called entanglement with vacuum, which 
turns out to be a notion representation dependent [49]. Number operator n^R' , fc') is 
added explicitly as an element of the Lie algebra as there exist representations where 
n b (R,k) / a b (R,k)^ a b (R,k). 
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Let us note that e~ lk ' x g a ® {R, k) is accompanied by the creation operator ao(R, ky . 
Had we employed in this place an annihilation operator, as it is typically done in the lit- 
erature, we would have to resign either from hermiticity of the 4-potential or positivity 
of the scalar product in the space of states. Both cases lead to nonunitary evolutions. 
Our formulation leads to a manifestly covariant quantization procedure with Hermi- 
tian 4-potential and does not involve unphysical indefinite-metric 'Hilbert space' [401 
IIP] , The field commutator 

[A a (x),A b (y)] = ig ab D(x - y) (206) 
involves a generalized Jordan-Pauli 'function' 

D(x) =i J dkdRI(R,k)(e~ tkx -e lhx ). (207) 

The fact that we work with twistor-like spin-frames guarantees that the four operators 
occurring in (1201[) - (|204p split into two different massless representations of inhomoge- 
neous SL(2,C): ai(R,k), a,2(R,k) (spin 1), and ao(R,k), a^(R,k) (spin 0). Transfor- 
mation x a >-> A a b x b + y a has to be represented unitarily at the operator level, i.e. there 
must exist U (A, y) such that 



U(A,y)^a(±,R,k)U(A,y) = e lp - y e ±2l0{A > k) a(±, A' 1 R, A _1 fc), (208) 

U(A,y^a 3 (R,k)U(A,y) = e lp y a^A' 1 R, A' 1 k) , (209) 

U(A, y) f a (-R, fc) f t/(4 y) = e ip ' v ^(A' 1 R, yl" 1 *;^. (210) 

Alternatively, 

U(A,y)U a (x)U(A,y)=A a b A b (A- 1 (x-y)). (211) 



For twistor-like spin-frames the above two sets of transformations are not inconsistent. 

Electromagnetic qubits are associated only with the spin-1 part (|208p . How to 
explicitly construct U (A, y) depends on the choice of representation of the oscilla- 
tor algebra. In general, formulas involving irreducible representations with I(R, k) = 
Z — const lead to mathematical inconsistencies typical of standard quantum field the- 
ory. 

In order to define a mathematically precise representation one has to carefully define 
Dirac delta, and take a nontrivial I(R,k). This type of quantization is nonstandard 
but I believe it is exactly what should be done (to understand how it works for fields 
interacting with sources, see [10|). 



6.2 Reducible representations of harmonic-oscillator Lie algebra and transformation 
properties of field operators 

One begins with four operators, oq, oi, a2, 03, satisfying commutation relations typ- 
ical of an irreducible representation of CCR: [a a ,ajj = <5 a t>l- Let j0) denote their 
common vacuum, i.e. a a |0) = 0. Now take kets |fe) and \R), normalized with respect 
to SL(2,C)-invariant M-shaped delta functions (see Appendix), (k\k'} = <5o(fc,fc')i 



32 



(R\ R ) = Si(R.R'). The reducible representation is parametrized by a natural num- 
ber N. For N = 1 the Hilbert space, denoted by is spanned by kets of the 
form 

fl,*,no,ni,na,n8> = \R, k) ® v °' \ 1 ' 2 -L_. 3 > |0), 

l-R, fe) = \R) ® |fe)- The 1-oscillator representation is defined by 

a a (R,k, 1) = \R,k)(R,k\ <g> a a , (212) 

n a (-R, fc, 1) = l-R, fc)(-R, fe| ® a-Ua, (213) 
I(R,k,l) = \R,k}(R,k\®l. (214) 

Operators /(ii, fc, 1) form a resolution of unity 

dkdRI(R,k,l)=I®I®l = I(l), (215) 

Let us note that (|213p is constructed independently of a a (R, fc, 1). The form n a (R, fe, 1) = 
a a (R, fc, l)^a a (i2, fc, 1) is possible if one works with M-shaped Dirac deltas 'regular at 
zero', i.e. 6a(k,k) = Sx(R,R) = 1 (see Appendix). Quantized analogs of the wave 
functions occurring in (|199|l . 



o a (fc,l)= / dRa a {R, fc, 1) = J® jfc)(fej ® a a , (216) 
are in this representation well defined. Now, let 



Pa(l) = 


—I (Si J dk k a \k) (k\ ® a a a a 




(217) 




J dkk a (n 1 (k, 1) + n 2 (k, 1) + n 3 (fc, 1) - 


- n (fc, 1)), 


(218) 




y dfc k a (n+ (fc, 1) + n_ (fc, 1) + n 3 (fc, 1) 


- «o(fc, 1)), 


(219) 


n a (fc, 1) = 


y dRn a (R, fe, 1) = 1 ® |fc)(fe| ® a a a a , 


a = 0, 1,2,3,+,- 


(220) 


a± = 






(221) 


h = 


i(a|a2 — a^ai) = a^a+ — a^a_. 




(222) 



The representation of inhomogeneous SL(2,C) then reads 

• ^-/dM^K^UTH,.*****. (223, 

(7(l,»,l)_e« ,p -<». (224) 

One immediately checks that 

U(A, y, l) f a(±, R, fc, 1)C/(A, y, 1) = e lp ' v e ±2l0[A ' k) a{±, A~ X R, A~ x k, 1), (225) 

U(A, y, l^a 3 (R, fc, 1)U(A, y, 1) = a^A' 1 R, A~ x k, 1), (226) 

[/(A, y, l) f a (R, fe, l) f ^(A y, 1) = e ip ' y ^(A' 1 R, A^k, l) f , (227) 

U(A, y, l)U a (x, 1)U(A, y, 1) = A a b A b (A~\x - y), l). (228) 



33 



Moreover, splitting the 4-momentum into electromagnetic and scalar parts, 

Pa(l) = P a {12} (l) + P Q {3} (1) - P a {0} (l), (229) 
Pi 12} W= Jdkk a (n+(k,l)+n-(k,l)), (230) 

Pi 3} W = J dkk a n 3 (k,l), (231) 

Pi° } (l) = J dkk a n (k,l), (232) 

we note that the three types of 4-momenta are not mixed with one another by SL(2,C) 
transformations: 

U(A,y,rfpV 2} (l)U(A,y,l) = j dfcfc a (n+(^- 1 fe, 1) +n-(A~ 1 k,l)) 

= A a b P b {12 \l), (233) 
[/(yl,y,l) t P a {3} (l)[/(yl,y,l) = Jdkk a n 3 (A- 1 k,l)=A a b Pl 3} (l), (234) 

U(A,y,l^P^ 0} (l)U(A,y,l) = J dk fc a n (yl- 1 fe, 1) = /l a fc P fc {0} (l). (235) 

For arbitrary N the representation is constructed as follows. Define 

H(N) =H(1)®...®H(1) = H(1)® N (236) 

s v ' 

N 

and let A be an arbitrary operator defined for N = 1. Let 

A {n) = 7(1) 7(1) ®A (8) 7(1) ® . . . ® 7(1) . (237) 

S v ' S v ' 

n— 1 TV— ri 

The N oscillator extension of a a (R, k, 1) is defined by 

N 

asL (R,k,N) = -j=^2a a {R,k,l) {n \ (238) 

n=l 

AT 



n a (R,k,N) = ^n a (7i,fc,l) ( " ) , (239) 

71=1 

AT 

7(7J,fc,iV) = 1^7(R,fe,l)<"), (240) 



N 

n=l 

U(A,y,N) = U(A,y,lf N . (241) 

Note that for N > 1 n a (R,k,N) a a (R, k, N)^a a (R, k, N). The 4-momentum de- 
scribes N noninteracting 4-dimensional oscillators of indefinite frequency 

N 



fa(A0 = $^Pa(l) (n) (242) 
n=l 

= Jdkk a (n + {k,N)+n-(k,N) + n 3 {k,N)-n {k,N)). (243) 
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The latter formula justifies the definition of n a (-R, k, N). I(R, k, N) again form a res- 
olution of unity 

J dkdRI(R,k,N) = I(N) = I{1)® N . (244) 

Of course, 

U(A, y, iV) f a(±, R, k, N)U{A, y, N) = e lp ' v e ±2l0(A ' k) a(±, A' 1 R, A~ x k, JV)(245) 
U(A, y, N^a 3 (R, k, N)U(A, y, N) = e lp ' y 'o 3 (A _1 R, A~ 1 k, N), (246) 
U(A,y,N^a (R,k,N^U(A,y,N) = e lp ' y aQ (A~ 1 R,A~ 1 k,N)' i , (247) 

and 

U(A,y,N^A a (x,N)U(A,y,N) = A a b A b (A-\x-y) J N), (248) 

so the representation of inhomogeneous SL(2,C) again splits into spin-1 and spin-0 
parts in the Fourier space, but at the level of Minkowski space A a (x, N) is an ordinary 
world-vector field. 



6.3 Polarization operators for fields quantized in reducible representations of 
harmonic-oscillator Lie algebras 

Our goal is to define polarization operators that have transformation properties anal- 
ogous to those based on Pauli-Lubanski vectors, but that are nevertheless linear in 
number-of-particles operators. In order to do so, let us note that for N = 1 there exists 
an analog of first-quantized 4-momentum, 



(249) 
(250) 



K a = I® / dkk a \k){k\ (g) 1, 
U(A, y, l)^K a U{A, y, 1) = A a b K b . 



This operator naturally occurs if one starts with 4-momentum that includes vacuum 
energy (see Appendix). Indeed, 

Pa{l)+K a =P (l)+/(g> Jdkk a \k){k\®^(l + l + l-l) 

= -J(g> Jdkk a \k){k\®^(aia a + a a aiy 

K a is in this representation well defined and commutes with elements of the oscillator 
Lie algebra. Generator 



M ab (l) = i 



dU(A, 0,1) 



dy ab 



y ab =Q 



= L ab (l) + S ab (l) 



(251) 



splits into orbital and spin parts, where 

SabW = I 



J dk\k)(k\®iJLe 2ie ( A M J * 



= 21® I dk\k){k\i-J?-rQ A ,{k)Air A ' (fc) 



Q y ab 

d 

,,ab 



y ab =0 

i 

y ab =Q 



■h 



(252) 
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and uiA(k) is any spin-partner of ir A (k). Now consider 

W a (l) = K b (l)*S ab {l) =7® [ dkk a \k)(k\®J 3 (253) 



dk 



k a (n+(k,l)-n-(k,l)\. (254) 



(255) 



An extension to A^ > 1, which is linear in n a (fc, N), satisfies 

N 

W a {N) = ^ W a {l) {n) = / dkk a (n+(k,N)-n-(k,Ny). 

n=l 

This type of extension is analogous to the relation between -Pa(l) and P a [N). 

6.4 Linear polarizations 
For N = 1 let us define 

V g (l)=I® I dk\k){k\ ® e 8(fe)(^«2-^ai) ( 25 g) 



satisfying 

Vg(l)oi(fe, l)V e (l) t = ai (fe, f) cos<9(fc) - o 2 (fe, f) sin6>(fc) = ag(k, 1), (257) 
V e (l)o 2 (fe, l)V e (l) t = o 2 (fc, f) cos<9(fc) + oi (fe, 1) sin0(fc) = o e /(fe, 1). (258) 
For arbitrary A^ > 1 we define 

V e (N) = Vg(lf N , (259) 
Ve(AT)ai(fc,iV)T/ e (iV) t = ai(k,JV)cos0(fc) - a 2 (fc, AT) sin 6>(fc) = ag(k,N), (260) 
Ve(Ar)a 2 (fe, A r )V e (A r ) t = a 2 (k,N) cos 9(k) + ai(fc, AT) sin 6>(fc) = a 9 / (fe, AT) ,(261) 

V 9 (A>(±,fe, A')'^(A') t = e ±W[k) a{±, fe, AT) = (a e (fc, AT) ± i a e , (fe, AT)) (262) 

Number operators corresponding to linear polarizations are defined analogously 

Vg{N)m(k,N)V s {N^ = n { (fe,JV), (263) 
Vg(N)n 2 (k,N)Vg(N)^ = ng,(k,N). (264) 

A yes-no observable associated with linear polarization can be defined as 

Y (h,N)=n e (k,N)-n ,(k,N). (265) 



7 Problem of relativistic analogues of EPR states of photons 

EPR states involve maximal entanglement in at least two different polarization bases 
(say, circular and linear). Such states can be easily invented also for photons, but 
the problem is with relativistic invariance of the EPR condition. At first glance the 
difficulty should not be a serious one since we know that the form (|28|l is SL(2,C) 
invariant: sab is a scalar representation of SL(2,C). However, for relativistic qubits 
the scalar-field condition has to be formulated at the level of unitary representation 
of SL(2,C), and it seems the problem is far from being fully systematized as yet. The 
basic question is if we can find scalar fields that involve maximal entanglement in two 
different polarization degrees of freedom, and in all reference frames. 
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7.1 Vacuum and multiphoton states 

To begin with, let us note that vacuum is any pure state annihilated by all annihilation 
and number operators. In reducible representations 

\0,N) = |O,1) 0JV , (266) 

|0,1) = J dkdRO(R,k)\R,k) ® |0), / dkdR\0{R,k)\ 2 = 1. (267) 

The analysis given in [10] suggests that the most natural choice of vacuum wave function 
is the separable state 0(R, fe) — Oo(k)Oi(R). Vacuum states are in this representation 
Bose-Einstein condensates at zero temperature. As such they are not unique, but the 
whole subspace of vacuum states is a relativistically invariant subspace of the Hilbert 
space of states. In field quantization based on twistor-like spin-frames an inhomoge- 
neous SL(2,C) transformation cannot produce particles, a fact which is not so obvious 
in general since Lorentz transformations in more standard approaches squeeze vacuum 
(a detailed discussion of this problem can be found in [10]). 

It is important to understand that the parameter N that characterizes a given 
reducible representation is unrelated to the number of photons. For example, 

n a {R,k,N)(a a (R,k,N)A J \0,N) = j (a a {R, k, TV) 1 ) 3 10, N) . (268) 

Weak (i.e. performed at the level of averages) limits TV — > oo reconstruct standard 
regularized formulas known from irreducible representations (for details see [40,10 

S3)- 



7.2 Scalar states 

The scalar-field condition means that 
>P(N) = J dRdR' dkdk'^(s,R,k;s',R',k')a(s,R,k,N^a{s',R',k',N^ 

(269) 

satisfies 

U(A, 0, N)*&(N)U{A, 0, TV) = $>(N). (270) 

Gauge independence implies that ^(s, R, fe; s' , R' , fe') can be constructed only by means 
of spinor fields tva (and thus does not depend on R and R') so that 

W = £ jaux s^'*(.M,: *>(.,*, 

. E /M* 1 ,»M,»Wtt'. (271) 

s,s — ± 

Since 

U(A, 0, /V) f a(s, fe, N^a(s', fe', N)^U(A, 0, N) 

= e -2is0(A,k) e -2is'0(A,k') a{s> A -l k Ar) t a(s ' j A -l k > t Ar) t ) 
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the problem reduces to finding li' satisfying 

9{a, k; a', k ') e -^{A,k) e -2is'0{A,k') = ^ A -i fc . ^ A - V). (27 2) 
Recalling that 

e - je (yl,fc)^ (fc) = ^^(A" 1 *), (273) 
e ie(yl ' fc) 7f^(fc) = A A , B '^(A^k) (274) 
we observe that only two types of contractions have the required form ((272]): 



(n A (k)7v A (k')ye- 2 ^ A ^e- 2ie ^ = (^(A^Tr^A"^')) (275) 
(7f A ,(fc)7f A '(fe')) 5 \^(^) e 2r0{A,k') = ^ A;( ^-l fe) ^' (A -l fe ')) 2 . (276 ) 
Since for two null vectors, k a ~ Tv A (k)n A i(k), l a = n A (l)n A i(l), we find 

fe- 1 = feX = 7r A (fc)7r A (i)7fA'(fe)7r A 'W = \^ A {k)-K A {I) | 2 , (277) 
the expressions (7^/ 

(fc)7f A '(fe')) and 7r yl (fc)7r A (fc') are proportional to each other, 
the proportionality factor k ■ k' being SL(2,C) invariant. 
So the corresponding solution reads 



<P(N) = J dkdk' F + (n A (k)-K A (k'))a{+,k,N^a(+,k',N^ 

+ / dkdk' F-(n A ,(k)n A ' ' {k'))a(-,k, N) 1 a(-,k' ,N)\ (278) 



where F± are any functions satisfying F±(e l ^z) = e 2z ^ F±(z). In reducible representa- 
tions the integrals in (I278|) are well defined. 

Another solution of (f2"72")l is <F(±, fe; ±, fe') = 0, 

<F{±, fe; =F, fe') = c ± S (k, fe'), (279) 

W(N) = c / dfca(+,fe,iV) t a(-,fe,iV) t . (280) 



States generated by both types of 'F(N) are maximally entangled in circular polariza- 
tions. Now take a closer look at linear-polarization entanglement of (|278|l . Inserting 

a(±,fe,iV) t = -j=e ±ie ^(a e {k,N)^ia e ,{k,N)^, (281) 

into (E7HJ, 



9{N) = j dkdk' \ (F + (. . . )e i[*(feW)] +F _(. . . )e -W»)+*(fe')j) 
x (ag(k, N^aeik', TV ) f - a e , (fe, TV ) W (fe', iV) f ) 

+ J dkdk' i (f+(. . _ . . )e -i[fl(*)+«(*')]) 

x (a fl , (fe, N)^ag(k', iV) f + ao(fc, iV) W (fe', iV) f ) , (282) 
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we obtain a condition for maximally entangled linear polarizations, 

F-(n A ,(k)* A '(k')) =±F + (^(fc)^(fc'))e 2i[e( ' c)+e(fc ' )1 . (283) 
Combining this with (|278[) . we arrive at 

&(N) = J dkdk' F + (n A (k)n A (k')) 

x (a(+, k, AT)ta(+, k' , JV)+ ± e ^[e(fc)+e(fe')] 

(-,fe,iV) t a(-,fe',7V) t ), 

(284) 

which does not satisfy (|272[) . State ty(N)\0, N) is scalar but does not satisfy linear- 
polarization EPR condition. It seems there is no difficulty with extending the above 
argument to general elliptic polarizations. 

The second solution (|280[) is more interesting. Indeed, (I280[) can be written as 

*{N) = § / dk(ag(k,N^ag(k,N^ +ae'(k,N^ag,(k,N)^, (285) 

which is simultaneously maximally entangled in all linear polarizations. Let us note 
that although Q285P involves photons of the same momenta, this does not exclude 
applications to experiments with pairs of detectors localized in space in arbitrary lo- 
cations. 



7.3 EPR states involving different momenta 

Let us again consider the operator (|27f [) but with the kernel satisfying 

<P(s, k; s, k') = -&(s, k; s, k') = -&(s, k! '■ s , k). (286) 
Denoting k; -,fc') = i>(k,k')/2 = -^(k!,k)/2, we find 

9(N) = J dkdk' ^{k,k')a{+,k,N)^a{-,k' ,N) ] (287) 

= {a5-M) I dkdk' i/j(k,k')b 1 (k,N)h 2 (k' ,N)^ 



where 

bi(fe,JV)t\ = /o^V 1 (a(+,k,N)i 
6 2 (fe,AT)t J \ 7 8j \a(-,k,N^J' 

and the matrix is unitary and independent of fc. In particular, if 

-y S ) ~ ^2 lie - * 9 -ie i6 



(288) 



(289) 



the state <f r (A r )|0, N) is maximally entangled in both linear and circular polarizations. 
One can replace the antisymmetric condition (|286[1 by a symmetric one, but then we 
arrive at the state of the form (128411, 
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Let us now consider the issue of SL(2,C) invariance of linear-polarization entangle- 
ment of (|287[) . The transformed operator 



U(A,0,N)^(N)U(A,0,N) 

= / dkdk' ${Ak, Ak')e- 2l0{AAk) a(+, k, N^e 2ie{A - Ak ' ] a{-, k' , iV) f (290) 



is maximally entangled in all linear polarizations defined by operators 

(V 4e(fc) a(+, k, /V ) f + e W{k) a{-, k, TV ) f ) , (291) 

-j= (e- ie(fe) a(+, fc, TV) 1 - e l6{k) a{-, k, 7V) f ) , (292) 

where 9(k) — 9 + 2<9(yl, Ak). In general, components corresponding to different wave 
vectors get rotated by different angles. 

Putting it differently, the antisymmetry of ip(k, fc') is not conserved by SL(2,C) 
transformations since 

^(Ak, Afc ' )e -2ie(A,Afe) e 2i©(A,Afe') (293) 

is not antisymmetric in fc and k' , so that the state is no longer maximally entangled 
in the original linear polarizations. 

The latter effect is essentially the massless version of the Peres-Scudo-Terno phe- 
nomenon. In the massless case the polarization observables are not defined with respect 
to projections of the Pauli-Lubanski vector, and thus the argument based on w-spinors 
cannot be directly employed. On the other hand, Wigner phases depend only on di- 
rections of fc, so all parallel wave vectors correspond to the same rotation angle. In 
order to maintain maximal entanglement in EPR-type experiments, one has to em- 
ploy momentum-dependent polarization operators that compensate the presence of 
the Wigner angle 20 (A, Ak). 



7.4 Normalization of 2-photon states in reducible representation 

Let us pause here for a moment to discuss the issue of normalization of 2-photon states. 
As an exercise consider the simple scalar case F+(z) = z 2 , F— (z) = 0: 



dk dk' dl dl' [tt a , (1)ty A ' " (k A {k)-K A (fc ') 

x (0, N\a{+, I, N)a(+, l' , N)a(+, k, 7V) t a(+, fc', 7V) f |0, N) 
dk dk' dl dl' (n A , (1)tt A ' (Z')) 2 (n A (k)TT A (fc')) 

x (<5 (fe, l')S (k', I) + S (k', l')S a (k, Z)) (0, N\I(k, N)I(k', N)\0, N) 

(l - jj) fdkdk' (k-k') 2 \0 (k)\ 2 \Oo(k')\ 2 . 



2 
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Let us note that in irreducible representations, where I(k) — Zl, Z — const, an 
analogous calculation would involve the divergent integral 

Invariance of the norm under inhomogeneous SL(2,C) is a consequence of the fact 
that the vacuum wave function is a scalar field: Oo(fc) — > Oq(A l k) is equivalent 
to \Q,N) — > U(A,y, N)\0, N). (Inclusion of vacuum 4-momentum K a would imply 
Oo(fc) -> e^O^A^k).) 

Performing an analogous calculation for (I287|l . we obtain 



{0,N\&(N)i&(N)\0,N) = (l - j^j J dkdk' \^{k,k')\ 2 \0{k)\ 2 \0{k')\ 2 . (295) 

In reducible representations the state H/(N)\0, N), generated by (|280|l . is normalizable 
as well. Since this exercise is also instructive, let us explicitly compute the norm for 
N = 1 in two alternative ways (put c = 1 for simplicity). The first approach explicitly 
employs the form of the vacuum state and normalization of kets to M-shaped deltas, 

<P(1)|0,1) = J dka(+,k, l) f a(-,fc, l)^ J dk' dR' O 1 (R / )O {k / )\R / ,k / , 0,0, 0,0) 

= J dkdR\R,k)(R,k\ ® a^al J dk' dR 1 O x (R!)Oo(k' )\r! ,k' ,0,0,0,0) 
dkdROi(R)O (k)\R, k) ® a\_(J_\Q) , 

so that (0, l|«^(l) t «^(l) jO, 1) = / dR\Oi{R)\ 2 J dk\O {k)\ 2 = 1. The second way is 
based on properties of the M-shaped deltas (see Appendix 110. 5|l , and the fact that 
annihilation operators annihilate vacuum: 

(0,1^(1)^(1)10,1) 

dk / dk'(0, l\a(+,k, l)a(-,fc, l)a(+, fe', l) f a(-, fe', l) f |0, 1) 



lim lim lim lim I dk I dk' 

n\— foo tl2— >oo 713— foo 714 — >-oo 



x (0, l|a(- k, 1, JL)a(+, ft, 1, £)a(+, fe', 1, ^) f a(- ft', 1, ^) f |0, 1) 
= lim lim [ dk [ dk'(0,l\a(-,k,l,MS (k,k')I(k,l)a(-,k',l,M^\0,l) 

Til ^OO 77,4— >• OO / 7 

dfc(0, l\I(k, I) 2 j0, 1) = (0, 1| / dk I(k, 1)|0, 1) = (0, 1|0, 1) = 1. 



One analogously checks that, for arbitrary N > 1 and c = 1, (1280 fl implies 

(0,iV|!?(JV)t*(JV)|0,JV> = + f dk\Oo(k)f. (296) 
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8 EPR averages for linearly polarized photons — the reducible 
representation approach 

Let us turn to the question of EPR averages for quantum electromagnetic fields quan- 
tized in reducible representations of harmonic oscillator Lie algebras. The issue is im- 
portant since all the experiments performed so far were analyzed in terms of irreducible 
representations, so one might be tempted to conclude that the standard approach to 
quantization is supported by EPR-type predictions. 

Let us consider the EPR state ([287]). Yes- no observables are defined by (|265l) . The 
EPR average 

(0, N\9(N)^Y p {l', N)Y a (l, N)<P{N)\0, N) 
{0,N\&(N)^(N)\0,N} 

8 Q {l\l)\O {l)\ 2 $ dkm,k)\ 2 \Oa{k)? - lvKM')| 2 |Oo(0| 2 |Oo(0| 2 
v H) J dkdk'\^(k,k')\ 2 \Oo(k)\ 2 \0 (k')\ 2 

(297) 

involves sharp wave vectors I, l' , and is independent of N if N > 1 (^(JV^O, N) = 
for N = 1). Since localization of detectors leads to momentum spread, I € SI, I & Q , 
say, yes-no operators integrated over both sets 



Y a (N) = / dlY a (l,N), (298) 
J a 

Y'piN) = [ dl'Y p (l',N), (299) 



can be used to compute more realistic cases. For disjoint detectors, Q n O' 
(0, N\¥{N)^(N)Y a {N)9(N)\Q, N) 



(0,N\&(N)i&(N)\0,N) 



- cos 2(a - P)p nxn ,, (300) 



where 



POxO' = PQ'xQ 

_ 2 Sg£ J n ,dl , \i>(l,l / )\ 2 \Oo(l)\ 2 \Oo(l')\ 2 
f R3 dk f R3 dk' |V(fe, fe')l 2 |Oo(fc)| 2 |0 (fc')| 5 

_ 2 Jgdl J n ,dl'm,l')\ 2 X (l) X (l') 

J R3 dk J R3 dk' \i>(k, fe')| 2 x(fc)x(fe') ' 



(301) 



The form (IM!) employs the cutoff function x( k ) = \0(k)\ 2 /Z, Z = max fc |C»(fc)| 2 , 
< xC 8 ) 5= 1| whose appearance is typical of predictions based on the reducible 
representation. 

For identical detectors, Q = i? , 

(0, N\W{N)^Y'JN)Y a {N)^{N)\Q, N) 

MW3 , (JV)|0)Ar) =coB2(a-/9)p„ x(JJ ,_ n) . (302) 

The average (1302|l vanishes if J? = J?' = It? . In the most general case, with arbitrary 
overlap S2q = J? n Q , O = Q\ U J?o, Si = J?o U Si[, one finds 

- cos2(q - P)(pn 1 xn , 1 +Pf2 1 xn +Pn xn , 1 - Pn x(R 3 -n )) ■ ( 303 ) 
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Let us note that Bell's inequality can be violated only if 

" V2 



Pa 1 xn' 1 + POi x n a + Pn xn[ - Pn x (r 3 - n a ) > > ( 304 ) 



For example, let 

V(fe, fe') = /(fe)ff(fc') - f(k')g(k), (305) 

where supports Xt = supp(/) and X g = supp(g) are disjoint (say, photons of two 
different colors are emitted, like in parametric down-conversion experiments). If Xf C 

o, x g c n', n n n' = <t>, then 

_ 2 f^rfZ Jjy dl'\f(l)g(l')\ 2 X (l)x(l') + 2 J n dl J n , dl' \f(l')g(l)\ 2 X (l) X (l') 
Pnxn ' J n dk J n , dk> |/(fe)ff(fc')l 2 x(fe)x(fe') + So, dk J Q dk' \f(k')g(k)\\(k')x(k) 

= Pn 1 y.n , 1 +P(iixQ +Pa xn , 1 -Pa x(R 3 -o ) = ^ ( 306 ) 

so that the Bell inequality can be violated. If Xf U X g C (2 = J?' 



/ dZ'|/(Z) 5 (O| 2 x(0x(O+ / / dl'\f(l')g(l)\ 2 X (l)x(l') = 0, 

and no violation is found. 

In irreducible representations, where I(k) = Zl, we get identical formulas but with 

_ f n dif n ,di'M,i')\ 2 
pn ^' 2 J R3 dkJ R3 dk'Wk,k^- (307) 

In practice, probabilities Q30ip and (|307p are not easy to distinguish from one another 
if x(k) is sufficiently flat in the optical-range set of wave vectors. The EPR averages 
thus do not provide us with any practical clue about the choice of representation of 
harmonic oscillator Lie algebras appropriate for field quantization. 

Relativistic properties of EPR averages cannot be directly inferred from 1)30011 but 
one has to take into account the remarks made after (|290[) and repeat the whole cal- 
culation. I will not pursue the matter further here. 



9 Remarks on two related issues 

One can find in the literature a discussion of 'entanglement with vacuum' (where a 
single-photon state is entangled) and violation of Bell's inequality in vacuum (where 
the vacuum state is entangled). In light of what I have written above on EPR states 
an existence of such phenomena may seem weird, but the problem is purely semantic 
— different authors have different things in mind when they speak of 'vacuum'. 
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9.1 Entanglement with vacuum 

The idea ol entanglement with vacuum comes from a specific representation of the Lie 
algebra 

[dm.aj,] = 5nml, (308) 

where m, n are integers, and 

On = ..-l®l®0®l®l... (309) 

with a, [a, a 1 ] = 1, located on an "nth position". a n acts in the non-separable Hilbert 
space spanned by infinite tensor products of number sates jn), a} a\n) — n\n). Infinite 
tensor products make the set of basis vectors uncountable ('Orlov states' [48] . such as 
. . . |0) ® |3) g> |1) ® 1 4) ® |1) ® 1 5) ® |9) ® . . ., with numbers taken from the digits of tt, 
are indexed by real numbers). 

The identity at the right side of (|308f) is given by the infinite tensor product of 
identities 

1 = ...1®1®1... (310) 

while the vacuum is represented by 

|0) = . . . |Q) ® |0) ® |0)... (311) 

A single photon state that exists in a superposition of two different modes, say, 

(aj, + 4„) |0) = . . . |0) ® |1) ® . . . ® |0) ® |0) . . . 

+ ...|Q)®|Q>®...®|0)®|1>... 

resembles, up to all the problems with infinite tensor products, the 2-particle entangled 
state 

|1)®|0) + |0)®|1). (312) 

This representation is often treated (especially in the quantum optics literature) as the 
representation of the harmonic oscillator Lie algebra typical of quantum fields. 

We have seen, however, that EPR correlations of photons can be computed at 
a much more general level and do not rely on specific tensor product structures of 
Hilbert spaces in question. The reducible representations we have worked with allow 
us to speak of all the field modes, but involve only iVth tensor powers. Attempts of 
interpreting, say, 

(a(+,fc,A0 f +a(-,fc,A0 t )|0,A0 (313) 

as an EPR state do not make much sense. Moreover, one does not have to resort to 
reducible representations to show that a single-photon 2-mode state may not always be 
interpretable in terms of entangled states. It is sufficient (see [35]) to take the original 
(1932) representation of the Fock space |50II51| . 



9.2 Nonlocal properties of vacuum states in algebraic quantum field theory 



Summers and Werner [4] showed that vacuum states, equipped with all their properties 
assumed in axiomatic quantum field theory, are enough to violate the Bell inequality. 
The paper is simultaneously the first published account of relativistic Bell theorem 
(see, however, [3]). The results from [4] were subsequently generalized in a number of 
works (cf. [52] for a recent review). 

In order to grasp the main idea let us consider a simpler but physically simi- 
lar problem suggested in [53]. Consider two systems, 0\, O2, equipped with certain 
algebras A(0\), -4(0 2 ). We assume commutativity [^i,A 2 ] = if A\ £ A{0\) 
and A2 £ .4(02). The algebraic structure can be represented in a Hilbert space 
Ti. = Hi ® W2 where representatives of A\ and A2 are of the form A <g> 1 and 1 <g> B, 
respectively. The crucial assumption about a vacuum state \ f2) is its cyclicity with 
respect to both A(0\) and A(02)- What it means is that acting on \ fl) with operators 
of either A(0\) or .4(02) one can generate any vector in H. 

To make our analysis as concrete and simple as possible, assume that Hi and M.2 
are 2-dimensional. If vacuum | Q) is cyclic in the above sense, then all basis vectors of 
Ti. can be written as {A® 1)|^}- Let |0j), |lj), span Hj. It is obvious that \Q) cannot 
be a product state. However, any state from the Bell basis 

|!Pi) = -^(|0i)®|l a )±|li)®|0 a >), (314) 
|*±) = -4f|0i)® |0 2 )±|li)®|l 2 >), (315) 



can play the role of \Q). Indeed, take A = |0i)(0i| - li)(li|, B = |0i)(li| + |li)(0i 
and \tt) = \<P+). Then 



(1®1)|«) = |^+), 

(A ® 1)|J?) = |*-_), 

(B ® l)\0) = |*+>, 

(AB ® = |$_). 



(316) 
(317) 
(318) 
(319) 



The same effect is obtained if one acts on the second qubit. Moreover, we can replace 
by any entangled state, say, a ground state of some 2-qubit Hamiltonian. \ f2) is 
then the lowest energy state which is cyclic with respect to commuting von Neumann 
algebras A(Oj). So this is precisely the vacuum in the sense of algebraic quantum field 
theory, but for a trivial toy model. The fact that 'vacuum' can maximally violate the 
Bell inequality is no longer weird. 

Our vacuum state |0, AT) is different. It belongs to a Poincare invariant subspace 
which is uniquely defined, but a single |0, AT) is neither unique nor relativistically 
invariant. |0, AT) is not cyclic either — it does not satisfy the axioms employed by 
Summers and Werner. In spite of that, the formalism leads to a well defined field 
theory, hopefully with no divergences, which is surprisingly close to the standard one 
in all the applications considered so far, simultaneously leading to new effects, testable 
at least in principle [101147] . 
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Fig. 1 The M-shaped function H320I I with a = 1, e = 1/2 (upper), and its Fourier transform 
(lower) . 

Appendices 

10 Dirac delta regular at zero 

10.1 M-shaped delta-sequences 

Let us consider the function shown in the upper part of Fig. 1. It is a particular 
example, for a — 1 and e = 1/2, of 




(320) 



(a > 0, e > 0). Its Fourier transform, 





(321) 
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lim S(x,a,e) = — , (322) 

is a real function shown in the lower part of Fig. 1 . The sequence S(k, a, i), with 
natural n (i.e. e = 1/n), is an example of what I call an M-shaped delta-sequence, and 
is in fact a particular example of a delta-sequence in the sense of [54]. Indeed, 

S(k, a, ±)dk = 1, (323) 

and for a function / left- and right-continuous at 



lim 

n— ¥00 



f(k)6(k, a ,k)dk= fi °- )+ 2 f ^ + \ (324) 

A peculiarity of M-shaped delta-sequences is their regularity at 0, 

S(0,a,i) = a, (325) 

for all n, so that 

lim 5(0, a i) = a. (326) 

71— >OC 

The fact that delta-sequences do not have to be divergent at the origin is, perhaps, not 
widely known but in itself is not new (an example of an analogous 'filtering function', 
vanishing at the origin, can be found in [55] . In what follows we restrict our analysis 
to M-shaped delta-sequences normalized by a — l/(2ir), 

5 M (k,±)=5(k,^,±). (327) 

Let us note that, in spite of regularity at 0, one finds 

roc 

lim / 8 M (k, ±)8 M (k, = co. (328) 



n—too 



The square of Dirac's delta thus will not exist even if we define delta in terms of 
M-shaped delta-sequences. However, 

lim lim / 5 M (k ±)5 M (k ±)dk = lim {„(0i)=i. (329) 



m— foo n— >-oo 



10.2 Plane waves and M-shaped delta-sequences 

We are heading towards an analysis of plane waves in terms of deltas that are regular 
at 0. To do so, we need delta-sequences that can be represented as scalar products of 
square-integrable functions. The simplest strategy is to start with the convolution of 
two M-shaped delta-sequences, 

/oo 
S M (k - k', ±)5 M (k>, ±)dk' = S* M (k, ±, kl (330) 
-oo 

poo 

jLirn^ / S M (k - k', ±)5 M (k', ±)dk' = 8 M (k, 1). (331) 
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The new sequence is again a delta-sequence, 



/oo />oo 
/ 6 M (k-k',±)6 M (k',±)dk>dk 
-CO J — OO 
fOO 

U 1 



M*'.^)d*' = l. (332) 



OO 



/OO 
/(*)ajf (*, £)d* 
-OO 

/CO /*CO 
/ f(k)SM(k-k',±)5 M (k',±)dk'dk 
-oo ./ — OO 

= lim / /(fc)fa(fc,l)dfc= /(0 - ) + /(0+) , (333) 

n— !-oo / Z 
«/ — OO 

but is not exactly M-shaped in the sense of the previous subsection. Indeed, 

/OO 
SM(0-k',±)S M (k',±)dk' 
-OO 

/CO 
<Mfc',i)<SM(fc',^)dfc' (334) 
-oo 

in general depends on n and m. The other properties are nevertheless analogous to 
M-shaped delta-sequences, 

1™ «M(0,i^) = M0,^) = l/27r, (335) 

m— >oo 

lim lim 5^(0, 1 i) = 1/2tt, (336) 



and 



Employing 



/OO 
«M(fc / ,£)fo(fc / ,£)<J* / = 00. 
-co 



1 f 00 

hh) = -k\ 5 *M(k, h iy kX dk = 2n6 M (x, ±)6 M (x, i) (337) 

J — OO 



we can write 



/OO 
-OO 



= 2tt/ 5 M (x,I)e lfea: «5 M ( a;; i)e" s3; dx 

«/ —CO 

= ^,^\k',±) = l-{k,±\k'^), (338) 
since 8m(x, jj) is real. Of course, (fc, ^} < oo for all fc, fc'. Now let us recall that 

lim 6 M {x,±) = ±, (339) 

n— >oo Z7T 
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and thus, formally, it is justified to write 

lim lim (*,I|*/,i>'=' JZ> e iik '- k)x dx, (340) 

while simultaneously we have shown that, for k = fc , 

lim lim (k, ijfc, -M = 2tt lim lim S* M (Q, ± -M = 1. (341) 



n— ^oo m— >oo n m n— >-oo m— foo 



There is no contradiction between the above two formulas — simply, integration does 
not commute with the limits n, m — > oo (integration must be performed first). 
A similar situation occurs with derivatives. Denoting 

(x|fe,i) = 2775m (x,kV kX (342) 

and taking into account 

Jim ii—fofol) = 0, AT = 1,2..., (343) 



we find 

lim (~) {x\k ±} = k N lim (x|fei) =k N e ikx . (344) 
Collecting all the formulas we have derived so far we are ready for generalization. 



10.3 Generalized function 8* M (k) 

From the point of view of the sequential approach to distributions [51] the sequences 
^A(k, i) = S(k,4n, i), 5 M (k, i) = 5(k, ±, i), and S* M (k, i, i) belong to the same 
equivalence class and thus define the same distribution. In consequence, any of these 
delta-sequences can be treated as a representative of the equivalence class and employed 
in computations involving Dirac deltas. 

However, we want to do more, and in particular want to include expressions such as 
Dirac delta at the origin. To do so, we will treat the limits lim n _>.oo lim m _>.oo 5* M (k, h,-k;) 
as a new type of generalized function, 8* M (k). The basic properties are as follows 



C 00 f{k)5* M {k)dk= lim lim /(fc)3 M (fc, ±)dk = J ^~ > | J{U+) ,(345) 

f^StrWdk = 1, (346) 

OmW 2 ^ = oo, (347) 

*Af(0) = ^, (348) 

5* M {x) = ± (349) 
Plane waves defined by 

(x\k) = lim (x|fc, i) = e ifer (350) 
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Fig. 2 The usual (/1-shaped) delta-sequences are special cases of M-shapcd delta-sequences 
— to see this one puts a = 4/e and defines 5^(k, e) = 5(k, 4/e, e). Here e = 1/2, a = 8 = 4/e 
(upper), and its Fourier transform (lower). 



satisfy 

(l±) N / x \k)= lim flA^/aHfc I) 
U dxJ v 1 7 n^oo V i dxJ x 1 n/ 



fe JV (a;|fe> ! (351) 



(fc|fc'> = lim lim <fc,i|fc',i) 

= 2tt5m(/c- k') (352) 
= lim lim {k,±\x){x\k',±)dx, 

(353) 

(fc|fc) = lim lim (k, hk, M = 2tt5* m (0) = 1 (354) 
n— >oo m— >oo 

These rules allow for all the standard computations involving the Dirac delta, but one 
can perform also certain new operations. For example, consider the expression 

c\k){k\k'){k'\y)f(k)dk, (355) 

where f(k) is, say, square-integrable, and which should be understood in the following 
sense 

lim lim lim lim [ (x\k, ±){k, ±\kf , ±){kf , ±\y)f{k)dk 
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7i2—^oo 713— too 



= Jim_lim_ / (x\k)2nS* M (k - k' , JL, ±)(k/\y)f(k)dk 



lim / (x\k)2n6 M (k-k',±)(k'\y)f(k)dk 



= 2-K{x\k'){k'\y)f{k') 

= 2n / (a;|fc)5(fc - k'){k'\y) f{k)dk (356) 



where the last formula was written in terms of the 'standard' Dirac delta, to make it 
more familiar, but we could replace S(k) by S^j(k). Simultaneously, 



ni— ^oo n-2— *-oc 77,3 — >oo 714— >-oo 



lim lim lim lim / ±){k, ±\y)f(k)dk 

= ( (x\k)(k\y)f{k)dk (357) 



The formulas reduce to 

\k)(k\k')(k'\ = 2ix5{k - k')\k){k\, (358) 

and 

\k){k\k){k\ = (359) 

that are mutually consistent if one treats each \k) as a limit of a separate sequence 
\k, i). Eqs. ()358p - (|359p can be conveniently written as 

\k){k\k')(k'\ = 27T(5M(fc - fc')|fc)(jfe|. (360) 



10.4 M-shaped deltas with respect to more general measures 

In this section we will concentrate on a generalization that is useful from the point 
of view of relativistic applications. Let us treat explicitly only a one-dimensional case. 
Let us assume that instead of dp we have to use a measure d/i(p) = p(p)dp, and an 
appropriate delta is needed, 

J dn(p')S„ M (p,p')f(p') = f( P ), (361) 

with 5um(p>p) — a i sa yi where a is a constant. The standard solution, 

8^P,p l ) = P {p'r 1 8{p-p l ), (362) 

if generalized to M-shaped deltas by 

SftM&p) = p(p')~ 1 Sm(p - p'), (363) 

implies <5^m(p>p) = p(p) _1 ^m(0) and will not lead to a independent of p. It follows that 
we have to proceed in a way different from (|363[) . So, let 8(p, a, i), 5(0, a, i) = a, be an 
arbitrary M-shaped delta-sequence discussed in the preceding sections. The sequence 



Sii(p,p',^) = p{p) 1 5(p~p',ap(p),±), 



(364) 
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p(p) = dfi(p)/dp, has the following properties 

lim / d»(p , )8 fl (p,p',±)f(j > ')= /Cp~) + /(p+) , (365) 

Wp>P> h) = a ( 366 ) 
and defines distribution <5 m m(PiP') satisfying 

drip'WM&pW) = + (367) 

<W(p>p) = a- (368) 



10.5 The meaning of products of operators occurring in harmonic-oscillator Lie 
algebra 

Consider N — 1 reducible representation of harmonic-oscillator Lie algebra. Define 

a b (R, k, 1, J_) = |H, J_ )(J?) J_| ® | fc) _L )(fc> J_| 06 . (369) 

The kets and bras are understood in the sense described in this Appendix. The products 
such as, say, a b (R,k, l)a c (R',k', l) 1 are understood in the sense of M-shaped Dirac 
deltas: 

a b (R,k, l)a c (R',fe',l) t 

= lim lim lim lim a b (R, k, 1, ±)oe{Bf, hf , 1, X, i)t. (370) 

tll->00 Tl2->OC tt',->00 fl^-KX) 1 2 1 2 

The order of limits is irrelevant. 



11 Noetherian construction of generators of N = 1 reducible 
representation — example of 4-momentum 

Consider any two fields A(x), B(x) satisfying d'Alembert equation. For example, let 
A(x) = A b (x, 1, = A c (x, 1, ^r), where the dependence on sequences 

means that the corresponding amplitude operators are constructed according to the 
recipe (|369|1 . d'Alembert equation implies conservation of the energy-momentum tensor 

T ab(x, & = -\{d a A{x)d b B{x) + d b A{x)d a B(x) - g ab d c A(x)8 c B (*)) , 

(371) 

i.e. d"T ab (x,±±lr,lr)=0, ! md 

P '(h J ^<*( x <>>*> (372) 
is independent of xq for all n\, . . . , n^- The limit 

Pa = lim lim lim lim P a ( — , — , -K- , 4-) 

= -7® Jdkk a \k)(k\0^(aia a + a a a^j = P a {l) + K a (373) 

is the 4-momentum, vacuum contribution included, corresponding to TV = 1 reducible 
representation of the harmonic oscillator Lie algebra. 
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